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Abstract 

Pilot-wave theories provide possible solutions to the measurement problem. In 
such theories, quantum systems are not only described by the state vector, but 
also by some additional variables. These additional variables, also called beables, 
can be particle positions, field configurations, strings, etc. In this paper we focus 
our attention on pilot-wave theories in which the additional variables are field 
configurations. The first such theory was proposed by Bohm for the free electro- 
magnetic field. Since Bohm, similar pilot-wave theories have been proposed for 
other quantum fields. The purpose of this paper is to present an overview and 
further development of these proposals. We discuss various bosonic quantum field 
theories such as the Schrodinger field, the free electromagnetic field, scalar quan- 
tum electrodynamics and the Abelian Higgs model. In particular, we compare 
the pilot-wave theories proposed by Bohm and by Valentini for the electromag- 
netic field, finding that they are equivalent. We further discuss the proposals for 
fermionic fields by Holland and Valentini. In the case of Holland's model we in- 
dicate that further work is required in order to show that the model is capable of 
reproducing the standard quantum predictions. We also consider a similar model, 
which does not seem to reproduce the standard quantum predictions. In the case 
of Valentini's model we point out a problem that seems hard to overcome. 
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1 Introduction 

In the pilot-wave theory of deBroglie and Bohm for non-relativistic quantum systems 
[1-3], also known as Bohmian mechanics [4], systems are described by both their wave 
function and by particle positions. The particles move along deterministic trajectories, 
under the influence of the wave function. The theory yields the same predictions as 
standard quantum theory, at least when the latter are unambiguous, when the particles 
are distributed according to quantum equilibrium.. The incorporation of particles in the 
description of quantum systems makes it possible to provide an objective description of 
the world, in which ambiguous notions such as "measurement" or "observer" play no 
fundamental role, so that there is no measurement problem. 

The particle positions are often referred to as "hidden variables" (like for example 
in Bohm's seminal papers [2,3]). However, as emphasized by Bell [5, p. 201], this 
terminology is rather unfortunate, for it is in the particle positions that one finds an 
image of the visible world. For example, whether Schrodinger's cat is alive or dead, 
or whether a certain outcome is observed in a measurement, depends on the way the 
particles are configured. The variable that is actually more hidden from us is the wave 
function, since it only manifests itself by its influence on the particles. 

Bell instead preferred to call the particle positions beables [5] . The beables represent 
the ontology of the theory, i.e. the things which are assumed to exist, independently 
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of being observed in any way. 3 The beables of a theory need of course not be particle 
positions, they could for example be fields or strings. 

When it comes to quantum field theory two possible choices of beables immediately 
come to mind, namely particle positions and fields. These choices are suggested respec- 
tively by the Fock representation and the functional Schrodinger representation (where 
states are functionals on a space of fields) of standard quantum field theory. In this 
paper we restrict our attention to fields as beables. The first pilot-wave theory with a 
field configuration as beable was given by Bohm already in 1952 [3]. In the appendix 
of his paper he considered this approach for the quantized electromagnetic field. Since 
Bohm, there have been several proposals for pilot- wave theories for quantum fields where 
the beables are fields. The purpose of this paper is to present an overview and further 
development of these proposals. Our attention is focused on the systematic development 
of the pilot-wave approach, rather than the study of it in various physical situations. 
For the latter we refer to [9-20] . 

The first quantum field theory we consider is the bosonically quantized Schrodinger 
field (which was discussed before in [17,21,22]). This is a simple field theory which 
serves as a first illustration of the techniques used in developing a pilot-wave theory. 

Then we consider the free quantized electromagnetic field. First, we show that 
Bohm's pilot-wave model can be derived by following a general recipe according to which 
beables are only introduced for gauge independent degrees of freedom. Bohm himself 
started with imposing the Coulomb gauge. There is also a different approach to gauge 
theories by Valentini (which was applied to the free electromagnetic field in [23,24], 
scalar quantum electrodynamics in [24] and the abelian Higgs model and non-Abelian 
Yang-Mills theories in [25]). Instead of first identifying the gauge independent degrees 
of freedom or imposing some gauge fixing, Valentini introduces also dynamics for some 
gauge degrees of freedom. The actual beable is then considered to be an equivalence 
class of fields related by the gauge transformations. We will compare the approach of 
Bohm and that of Valentini in detail for the case of the free electromagnetic field. We 
find that while the ontologies are strictly speaking different, the theories are empirically 
indistinguishable, even in quantum non-equilibrium. 

After discussing the pilot-wave approach for the free quantized electromagnetic field, 
we present a pilot- wave model for scalar quantum electrodynamics (a scalar field coupled 
to the electromagnetic field) and the Abelian Higgs model. In each case we only intro- 
duce beables for gauge independent degrees of freedom, thereby avoiding the surplus 
structure in Valentini's approach. The Abelian Higgs model provides a simple example 
of a theory that exhibits spontaneous symmetry breaking and the Higgs mechanism. 
We will see that the pilot-wave approach yields an elegant description of this. 

We also briefly consider non-Abelian Yang-Mills theories, which are the natural gen- 
eralizations of the electromagnetic field theory. However, due to technical complications 

3 More precisely the particle positions are local beables [6, 7], because they relate to particular regions 
in space-time. In the case of a field ontology the local beables are the fields restricted to space-time 
regions. The wave function is also a beable but a non-local one. In the terminology of Diirr et al. the 
particles constitute the primitive ontology [4, 8]. This notion is closely related to that of local beables. 
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we do not provide a pilot-wave model for these theories. The technical complications 
have to do with the difficulties in identifying and isolating the gauge independent degrees 
of freedom. We expand on this in a short note. 

Then we turn to fermionic quantum field theories. There are two suggestions for 
introducing field beables for fermionic quantum fields. There is Holland's suggestion 
[17, 22] to take a field of angular variables as beable and Valentini's suggestion [23, 24] 
to introduce anti-commuting fields as beables. We discuss in detail the ontology of 
Holland's model and point out that it is not clear whether the model is capable of 
reproducing the standard quantum predictions. We also consider an alternative model, 
which can be developed using the same techniques as in Holland's model. For this model 
it is simpler to analyse whether or not it reproduces the standard quantum predictions. 
We find that it seems incapable of doing so. In the case of Valentini's model, we point 
out a problem that seems hard to overcome. We find that although it is possible to 
introduce a dynamics for anti-commuting fields, it is unclear what the actual guidance 
equations and the associated equilibrium probability distribution should be, in order to 
have a pilot-wave model that reproduces the standard quantum predictions. 

Finally, we briefly discuss the pilot-wave models proposed by Struyve and Westman 
for quantum electrodynamics [26,27]. The first model is radically minimalist because it 
contains only beables corresponding to bosonic degrees of freedom, the electromagnetic 
field degrees of freedom in this case, and none for the fermionic degrees of freedom. 
The second model is an extension of the first with additional beables for the fermionic 
degrees of freedom, such as for example a field beable corresponding to mass or charge 
density. 

In this paper we restrict our attention to field beables. However, we want to empha- 
size that a field ontology is not the only possible pilot-wave ontology for quantum field 
theory. One possible alternative is a particle ontology, which seems to be better suited 
for fermionic field theories. There is for example the seminal work of Bell [28] for lattice 
quantum field theory where the beables are the numbers of fermions at each lattice 
point. This work has recently been generalized in two different ways to quantum fields 
on the continuum. One possible continuum generalization was presented by Diirr et al. 
[29-33]. In this approach particles travel deterministically along trajecties, interrupted 
by stochastic jumps which are usually associated to the creation and annihilation of 
particles. A different continuum generalization was presented by Colin [34-36] for quan- 
tum electrodynamics, and generalised by Colin and Struyve [37] for other interactions 
in the standard model. In this approach the Dirac sea is taken seriously, in the sense 
that particle positions are introduced that correspond to the Dirac sea. A difference 
between the two models is that the model of Colin is deterministic, whereas the model 
of Diirr et al. is, just as Bell's lattice model, stochastic. The fact that there are two 
generalizations of Bell's lattice model for the continuum originates in a different reading 
of Bell's work, see e.g. [37,38]. 

Before discussing the various pilot-wave models for quantum fields, we will first recall 
the pilot- wave theory of deBroglie and Bohm for non-relativistic quantum systems. 
Various aspects of this pilot-wave theory, like the explanation of how the standard 
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quantum predictions are reproduced, are straightforwardly generalized to the case of 
field beables. 

We also outline in general how field beables are introduced. We thereby start with 
recalling how a quantum field theory can be obtained by quantizing a classical field 
theory. There are different ways of quantizing a theory, but we focus on a way that 
is particularly convenient for developing a pilot-wave theory. In particular, it is inter- 
esting to see how gauge theories are quantized, since this will also reveal the status of 
gauge invariance in the corresponding pilot-wave models. For clarity, these quantization 
methods are standard in quantum field theory and have a priori nothing to do with 
pilot-wave theory. 

Throughout the paper we use natural units in which h — c — 1. 

2 From a classical theory to a pilot- wave theory 

There exist several procedures to construct a quantum theory starting from a classical 
theory, see e.g. [39] for a recent overview. These procedures are called quantization pro- 
cedures. Some well-known examples are canonical quantization, path integral quantiza- 
tion and BRST quantization. Even within these different procedures, one has different 
approaches. The different quantization procedures usually yield, though not always, 
equivalent quantum theories. 

For the construction of a pilot-wave model for the bosonic quantum field theories we 
consider in this paper, we found it convenient to start with a formulation of the quantum 
theory in the context of a particular approach to canonical quantization. This approach 
consists in identifying the "unconstrained degrees of freedom" and by quantizing only 
these degrees of freedom. Below we make clear what exactly is meant by this. Just to 
give an idea we can already mention that in the case of gauge theories the unconstrained 
degrees of freedom are simply the gauge independent degrees of freedom. Beables are 
then introduced only for some of the unconstrained degrees of freedom. 

For fermionic fields the canonical quantization can be carried out in a similar way, 
but the step of introducing beables causes complications. In this section we will therefore 
only consider bosonic field theories. The discussion of fermionic field theories is deferred 
to section 9. 

Although our approach works well for the field theories we consider in this paper, it 
is not necessarily applicable to other field theories. In such cases it could be interesting 
to look at alternative quantization procedures. Different quantization procedures may 
further lead to different formulations of a quantum theory, which in turn may suggest 
the construction of different pilot-wave models. 

In order to illustrate the procedure of canonical quantization, without the techni- 
calities and subtleties that go with quantum field theories, we first consider theories 
with a finite number of degrees of freedom. After that, we discuss how the transition is 
made to a field theory. We start with explaining the particular approach to canonical 
quantization. Then we consider how to construct a pilot-wave model for the obtained 
quantum theories. We also give a short discussion of the pilot-wave theory originally 
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presented by de Broglie and Bohm for non-relativistic quantum systems. 

2.1 Finite number of degrees of freedom 
2.1.1 Canonical quantization 

We first consider theories that concern a finite number of degrees of freedom. We assume 
that they can be described by the Lagrangian formalism. The Lagrangian L(q, q) is 
then a function of some coordinates q n , n = 1, . . . , N, and the corresponding velocities 
q n = dq n /dt. The corresponding Euler-Lagrange equations of motion read 

^-^ = 0. (1) 
dtdq n dq n 

In the canonical quantization approach, the first step is to make the transition from 
the Lagrangian formulation, where the dynamical variables are the velocity phase-space 
variables q n and q n , to the Hamiltonian formulation, where the dynamical variables 
are the momentum phase-space variables q n and p n . Canonical quantization then pro- 
ceeds by associating operators to the momentum phase-space variables and by imposing 
certain commutation relations for these operators. 

In order to arrive at the Hamiltonian formulation we first need the momenta p n that 
are canonically conjugate to the coordinates q n . These are defined as 

dL 

Pn = W M- (2) 

In case these relations can not be inverted to yield the velocities q in terms of the 
coordinates q and the momenta p, the Lagrangian is called singular. The Hamiltonian 
picture for singular Lagrangians requires some care. The formalism that deals with such 
Lagrangians was originated by Dirac [40]. Extensive reviews on Dirac's formalism can 
be found in [40-44]. Here, we just give some results that are relevant for the paper. In 
order to keep the presentation simple, we ignore subtleties that are not important for 
the theories that we will consider here. 

The fact that the relations (2) are not invertible implies that there are certain re- 
lations between the momentum phase-space coordinates (q,p). These relations can be 
represented as Xm(q,p) — 0, m — 1,...,M', and are called constraints. There may 
be further constraints that arise from the consistency requirement that the constraints 
Xm(<2SP) = are preserved in time. There is a well-defined algorithm to find these extra 
constraints, but we do not present it here since the constraints are well known for the 
theories we will consider. We just assume that we have a complete set of independent 
constraints Xm{c[,p) — 0, m — 1, . . . , M (where M' ?C M). 

One can further define the Hamiltonian function 

H = ^ VrAn - L , (3) 

n 
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which is a function of the coordinates q and the velocities q. By making use of the 
relations for the momenta (2) and the equations of motion (1), it can be shown that 
the Hamiltonian can be written as a function H(q,p) of the momentum phase-space 
coordinates (q,p), despite the fact that some of the velocities q n are not expressible in 
terms of the coordinates and the momenta. 

For singular systems the equation of motion for a dynamical variable F(q,p) is not 
given by the usual relation 

F=[F,H] P , (4) 

where [., ]p is the Poisson bracket. The reason is that the constraints have to be taken 
into account. There exists a generalization of the above equation of motion, but for the 
particular quantization procedure we will consider here, we do not need to go into this. 
A constraint Xm is called first class if 

[Xm,Xm']p =0, m' = l,...,M, (5) 

Xl=-=XM=0 

otherwise it is called second class. Theories with first class constraints and theories with 
second class constraints are quantized differently. While one can have theories with both 
first and second class constraints, we only consider theories which have either of those. 



First class constraints If the constraints are first class then we are dealing with a 
gauge theory. A gauge theory is hereby understood as a theory in which the dynamics 
of some variables is not determined by the equations of motion and the specification of 
initial data. 4 In a theory with first class constraints this underdetermination arises as 
follows. The equation of motion for a dynamical variable F is given by 

M 

F=[F,H] P + J2[F,Xm}pu m , (6) 

m=l 

where the u m are arbitrary phase-space functions which may depend on time, and where 
the constraints may only be taken into account after the Poisson brackets have been 
evaluated. If the dynamics of some dynamical variable depends on those arbitrary 
functions u m , then its time evolution is underdetermined, since different choices of the 
u m will lead to different evolutions in time. A map that maps solutions to solutions 
with the same initial data, as well as a combination of such maps, is called a gauge 
transformation. Solutions that are connected by a gauge transformation are called 
gauge equivalent. Dynamical variables for which the evolution does not depend on the 
arbitrary functions are said to be gauge invariant. In other words, for gauge invariant 
variables the time evolution is uniquely determined by the equations of motion and the 
initial data. These variables can be regarded as the physical degrees of freedom of the 
theory. From the equation of motion (6) it is clear that gauge independent variables are 

4 Although one can adopt alternative notions of gauge theories, this notion is very natural and 
commonly used for theories whose the equations of motion can be derived from a Lagrangian. This 
notion of gauge theories further plays an important role in standard methods of quantizing such theories. 
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those for which the Poisson brackets with the constraints are a linear combination of 
the constraints. 

The above definitions also apply in the Lagrangian picture. However, the gauge 
freedom shows up more explicitly in the Hamiltonian picture. 

Here we consider two ways to quantize a theory with first class constraints: 

• Quantizing the gauge independent degrees of freedom: It can be shown 
that a canonical transformation can be performed, at least locally, such that the 
new canonical variables can be written in terms of two sets q' and q", and their 
respective conjugate momenta p' and p", such that, in terms of the new variables, 
the constraints read p" = 0, see e.g. [43, pp. 36-45] or [45]. As such the variables 
q' and p' can be identified as the gauge independent variables, while the variables 
q" are gauge variables. 

The equation of motion for a function F(q',p') is generated in the usual way: 

F=[F,H ph ] P , (7) 

where 

HM,P') = H(q',p',q",p")\ pll=0 (8) 

and where H(q',p',q",p") is the Hamiltonian that is obtained from H(q,p) by 
the canonical transformation. The Hamiltonian H p h, which does not depend on 
the variables q" , is called the physical Hamiltonian. Subsequently, we will gener- 
ally omit the subscript "ph" of the Hamiltonian. The variables q" are the gauge 
variables; their motion is completely undetermined. Therefore the variables q", 
together with the variables p" (which are constrained to be zero), can be ignored 
in the description of the system. 

The theory can now be quantized just as in the case of a non-singular Lagrangian. 
Modulo some potential problems, which we touch upon below, this can be done as 
follows. First one associates operators <f and p' to the canonical variables q' and 
p', and imposes the commutation relations 

[Z,%\=S5 n . (9) 

Then, using the Schrodinger representation 

d 

q r ^ q r , p r ^ -1-^7- , (10) 

where the operators act on square integrable wave functions ip(q f ), for which the 
inner product is given by 

<lMlfe> = / dq'rMMl') , (11) 
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one obtains the Schrodinger equation 



1 dt 



H^q',t), 



(12) 



where H is the operator associated to H ph . 

This quantization procedure appears very natural, but there might be difficulties 
in applying it. Although these difficulties do not appear in the theories we consider 
here, except maybe in non-Abelian Yang-Mills theories, we want to briefly mention 
some of these here. For more complete discussions see [41-44]. 

The first difficulty that may arise concerns the identification of the gauge inde- 
pendent degrees of freedom. The problem of finding a canonical transformation 
which separates the gauge degrees of freedom from the physical degrees of freedom 
has so far not been constructively solved (except in some simple cases). It is also 
possible that the canonical transformation can not be performed globally. But for 
the other theories we consider here, it will always be easy to find a global canonical 
transformation. 

Secondly, even if the gauge independent degrees of freedom can be identified, it 
might be that their phase-space is such that the above quantization procedure can 
not be applied. When the phase-space of gauge independent degrees of freedom, 
which is parametrized by q' and p', is given by M. 2 ( N ~ M ^ there is no problem, but 
in the case of other phase-spaces it might be necessary to modify the quantization 
procedure [46] (this is for example the case for a particle moving on the positive 
half-line, where one needs to impose different commutation relations). However, 
for the theories we consider in this paper we do not need to worry about this; 
modulo the problems that appear when dealing with an infinite number of degrees 
of freedom, we can always quantize as outlined above. 

It can further be shown that the physical degrees of freedom q' and p' are unique 
up to a canonical transformation [43, p. 40]. Choosing to quantize a different set of 
canonical coordinates which parametrize the physical space does not always lead 
to an equivalent quantum theory. Possible inequivalence might arise due to the 
operator ordering ambiguity, i.e. equivalent orderings of the classical coordinates 
might yield inequivalent orderings of the corresponding operators. However, such 
ambiguities do not concern us here. 

In presenting pilot-wave models for gauge theories, we will always start with pre- 
senting the quantum field theory in terms of the gauge independent variables, 
according to the above framework. 

• Imposing constraints as conditions on states: A second way to quantize is 
to use the quantization rules as in the case of a non-singular Lagrangian 



[q r ,p s ] = i£, 



rs 



(13) 
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and to impose the constraints as conditions on states 



Xm\il>) = 0, m = l,...,M. (14) 

The states that satisfy the above condition are called gauge independent or physical 
states. 

This quantization procedure works when the commutators [£ m ', Xm"\ an d [xw> H] 
are linear combinations of the constraint operators £m (although such relations 
hold classically, quantum corrections might appear due to quantization [42,44]). 
The first condition is required for consistency, the second guarantees that physical 
states remain physical under the Schrodinger evolution. 

A potential difficulty with this approach lies in the definition of the inner product. 
Working in the Schrodinger representation q r — > q r , p r — > —id/dq r , the inner 
product 

(i/>iM = J dq^MMo) (15) 

will not be finite for physical states. In order to obtain a finite inner product the 
measure in the integral could be changed, for example by applying the Faddeev- 
Popov procedure, see e.g. [44, pp. 281-283]. (A simple example is given by the case 
where one of the constraints reads p r \ip) = 0. In the Schrodinger representation 
this implies that physical states tp(q) should not depend on q r , so that the above 
inner product becomes infinite. However, by inserting the distribution S(q r ) in the 
integral one obtains a finite inner product.) We will discuss this in more detail 
when considering this quantization procedure for the electromagnetic field. 

We mention this way of quantizing, because Valentini's pilot-wave approach fits 
in this framework. 



Second class constraints Just as in the case of first class constraints, the uncon- 
strained degrees of freedom can be isolated by performing a canonical transformation. 
It can namely be shown that a canonical transformation can always be performed, at 
least locally, such that the new canonical variables can be written in terms of two sets q' 
and q", and their respective conjugate momenta p' and p", such that in terms of the new 
variables, the constraints read q" = p" = 0, see [43, p. 27-35] and [45]. The equation of 
motion for a function F(q',p') of the unconstrained variables reads 

F=[F,H ph ] P , (16) 

where the physical Hamiltonian is given by 

H ph (q',p') = H(q',p',q",p")\ qf/=pll=0 (17) 

and where H(q',p', q",p") is the Hamiltonian that is obtained from H(q,p) by the canon- 
ical transformation. 
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The unconstrained degrees of freedom q' and p' are unique up to a canonical trans- 
formation [43, p. 32]. 

The resulting theory for the unconstrained canonical variables q' and p' can then be 
quantized in the usual way. 

Some comments on other types of quantization In the case of second class 
constraints the Dirac bracket can be introduced [40-44]. Unlike the Poisson bracket, the 
Dirac bracket is consistent with the constraints. This means that it makes no difference 
whether the constraints are imposed before or after the Dirac bracket is evaluated. 
An alternative way to quantize then consists in replacing the momentum phase-space 
coordinates by operators and by imposing commutation relations that are based on the 
Dirac bracket rather then on the Poisson bracket. However, we do not know how the 
resulting formulation could be useful for developing a pilot-wave model. Of course one 
could use the fact that the Dirac bracket can locally be expressed as a Poisson bracket 
for unconstrained variables, see [43, p. 30] and [45], but this then essentially amounts 
to the quantization procedure for systems with second class constraints we presented 
above. 

In the case of first class constraints one could also introduce extra constraints, also 
called gauge constraints, so that the total set of constraints becomes second class. One 
could then use techniques employed in the case of second class constraints to quantize. 
However, again we do not know how to exploit this feature for the purpose of construct- 
ing a pilot-wave model. In particular, it can easily be seen that a reformulation of the 
theory in terms of unconstrained degrees of freedom after gauge fixing, is equivalent to 
a such reformulation without gauge fixing. 

2.1.2 Construction of a pilot-wave model 

Consider now the construction of a pilot-wave model. For the moment we still consider 
systems that can be described by a finite number of degrees of freedom. The transition 
to a field theory is explained in section 2.2. 

Consider a classical Hamiltonian of the form 



where h rs is symmetric and where the phase-space coordinates (q, p) G M, 2N are uncon- 
strained (see [47] for the treatment of other Hamiltonians). The classical field Hamilto- 
nians that we will consider in this paper will have a similar form. Quantization proceeds 
by associating operators q and p to the canonical variables and by imposing the com- 
mutation relations 





r,s 



[q r ,Ps\ = i5 rs . 



(19) 



Using the representation 




(20) 
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the following Schrodinger equation is obtained 



. dij>(q,t) 
1 dt 



(4 £^<< + v <«>) *<«•'>■ (21) 



The wave function ip(q,t) can be regarded as the expansion coefficient (q\ip(t)) of the 
state \ip(t)) in the basis of states \q), where the states \q) are eigenstates of the operators 
q r , i.e. q r \q) = q r \q). 

There is an operator ordering ambiguity in the kinetic term (for example, regarding 
h rs as a metric, the Laplace-Beltrami operator could be introduced instead). However, 
the operator ordering that is used here is the most natural for the theories we will 
consider. With this operator ordering, the Hamiltonian operator is also Hermitian with 
respect to the inner product 

(H^) = J dqrftMMQ) ■ (22) 

Other operator orderings could potentially lead to different quantum theories. 

According to standard quantum theory, \ip{q,t)\ 2 dq = \(q\i[)(t))\ 2 dq yields the prob- 
ability that a measurement of the observable q yields an outcome which is in a volume 
element dq around the configuration q. The Schrodinger equation implies the following 
continuity equation for the density \ip(q, t)\ 2 : 

w + £f = °. (23) 



where 



dt ^ dq r 



is the probability current, with ip = l^e 15 . 

A pilot-wave model can be constructed by introducing a beable configuration q, 
which takes values in the configuration space M. N , and which evolves in time according 
to the guidance equation 

^ = ^ = E^- ( 25 ) 

This dynamics implies that if the beable configurations are distributed according to 
\ip(q, to) 1 2 over an ensemble at a certain time to, then they are distributed according to 
\ip(q, t)\ 2 at other times t. This means that the distribution \ip\ 2 keeps its functional form 
of the wave function tp. This property is called equivariance [4]. The distribution |-0| 2 
plays the role of an equilibrium distribution [4, 16, 23,48-50], similar to that of thermal 
equilibrium in classical statistical mechanics, and is called the quantum equilibrium 
distribution. In this paper we always assume quantum equilibrium. 
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In some models some additional beables will be introduced by means of Holland's 
local expectation value [17]. The local expectation value corresponding to a Hermitian 
operator Q and a wave function ip(q,t) is the function 

Qiq > t] - Re K%^yp ■ ( } 

A new beable Q can then be introduced, by evaluating the local expectation value for 
the actual value of the beable configuration q, i.e. 

Q(t) = Q(q(t),t). (27) 

As such its time evolution is completely determined by the time evolution of the wave 
function and that of the beable configuration. In quantum equilibrium, the distribution 
of the beable Q is given by 

p(Q,t) = J dq\^(q,t)\ 2 5(Q - Q(q,t)) . (28) 

This distribution is in general different from the distribution KQl^)] 2 , where the states 
\Q) are the eigenstates of Q. The definition by means of Holland's local expectation 
value only guarantees that the expectation value of the beable Q agrees with the quan- 
tum mechanical expectation value (ip\Q\ip). However, during actual measurements it 
might be that the beable Q becomes distributed according the quantum mechanical dis- 
tribution. This happens for example for such beables introduced in the context of the 
non-relativistic pilot-wave theory of deBroglie and Bohm, see [17, pp. 339-347]. Note 
that we definitely do not want to introduce a beable for every Hermitian operator in 
this way. In most cases the introduction of such additional beables is even unnecessary. 

We want to end this section with a word of caution concerning the choice of beable. 
The beables q r (and Q) that we introduced correspond to quantum operators q r (and 
Q) and one could get the impression that we could have chosen any operator for which 
to introduce beables. However, this is not correct. The beables should yield an image of 
the visible world (i.e. they should account for things like tables, chairs, cats, instrument 
pointers, etc.). Therefore, the beables should in the first place relate to an ontology 
in space-time. If the beables are for example associated to momentum operators or 
spin operators, without a specification of how these beables relate to, or influence the 
behaviour of, some ontology in space-time, it is unclear how they could yield an image 
of the visible world [51]. Here we are concerned with models in which the beables are 
fields in space and these could potentially yield an image of the visible world. 

2.1.3 Example: Pilot-wave theory of de Broglie and Bohm for non-relativistic 
quantum systems 

In this pilot- wave theory [1-3], the Schrodinger equation is given by 
.<9V>(xi,...,x n ,t) ( V? \ 

i Qj. = I - 2^ ^ + V(xi, • • • , x n ) I V(xi, . . . , x n , t) , (29) 
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where iftfai, ■ ■ ■ , x„, t) is, for all times t, a function on configuration space lR 3ra . The 
beables are n particle positions in physical space M 3 , denoted by x fc , k — 1, . . . , n, for 
which the evolution equation is given by 



1 



(^*v fc v-^v fc v*) 



(30) 



Xfc = 



2im k \ip\ 2 



The quantum equilibrium distribution is given by ■ ■ ■ , x„,t)| 2 . 

Bohm actually presented the dynamics in a Newtonian form, with an extra in- 
dependent potential, called the quantum potential, thereby regarding the guidance equa- 
tion as an extra constraint on the possible momenta. However, while this Newtonian 
formulation has some applications, it tends to obscure and complicate things. There- 
fore, we will consider only formulations in terms of guidance equations throughout this 
paper. 

Having the quantum equilibrium distribution for ensembles is a key ingredient in 
showing that the pilot-wave theory reproduces the predictions of standard quantum 
theory. As we will explain below, it guarantees that the Born rule is recovered. As shown 
by Diirr et al. [4], most initial configurations of the universe (relative to the natural 
measure |\l/(X)| 2 <iX, with \& the wave function of the universe) yield the quantum 
equilibrium distribution |?/>(x)| 2 for actual ensembles described by the effective wave 
function ip(x). As such, one also expects that most non-equilibrium distributions tend 
to evolve to equilibrium. This was indeed illustrated by numerical simulations [52,53]. 

Another key ingredient of the pilot-wave theory is that wave functions representing 
macroscopically distinct states (like for example a cat in the alive state or dead state, 
or macroscopically different orientations of an instrument needle) have their support 
approximately concentrated in different regions in configuration space. That is, these 
wave functions have their support approximately concentrated on configurations that 
on the macroscopical level yield the same familiar image of the macroscopic objects. 
As a result, in quantum equilibrium, the particles will typically attain one of those 
configurations and as such yield an image of the object. Even if the wave function 
happens to be in a superposition of macroscopically distinct states of a macroscopic 
object, the configurations will typically end up displaying either state of the object. For 
example in the case of Schrodinger's cat, where the wave function is in a superposition 
of the state of a live cat and a dead cat, the configurations will typically end up either 
displaying a live cat or a dead cat. In particular, in a measurement situation, the non- 
overlap of wave functions representing macroscopically distinct states of macroscopic 
pointers, like for example instrument needles, will guarantee that measurement results 
get recorded and displayed in the position configurations of those pointers. In addition, 
it will also give rise to an effective collapse. 

Let us now consider the pilot-wave description of a measurement situation in detail 
and see how the standard quantum predictions are reproduced. This analysis was first 
given in [3] and has been repeated many times afterwards, see e.g. [28,54] (for an 
extensive analysis see [55]). In standard quantum theory a measurement situation can 
be described as follows. Before the measurement the wave function is given by a product 



15 



i/j s i/j a , where ip s = X)j c »V'i * s the wave function of the system, with ^| the different 
possible eigenstates of the operator that is being measured, and ip a is the wave function of 
the macroscopic pointer in the ready state (one could also include the rest of the relevant 
environment in if) a ). The macroscopic pointer could for example be an instrument 
needle. During the measurement process, the wave function evolves to the entangled 
form Y2i Ci^t^t- The ijjf are the different states of the macroscopic pointer corresponding 
to the different possible outcomes of the measurement (for example an instrument needle 
pointing in different directions). The wave function then collapses to one of the terms 
ipfipf, say i> s k il>l-, and this with probability \cu\ 2 . The wave function of the system ends 
up in the /c-th eigenstate of the operator and the wave function of the pointer represents 
the /c-th outcome of the measurement. 

In the pilot-wave theory the measurement situation is described as follows. As in 
standard quantum theory the wave function evolves into a superposition ^\ c^ipf, but 
this time there is no subsequent collapse. The measurement result does not get recorded 
in the wave function, but in the configuration (x", . . . , x°) of particle positions of the 
pointer. This can be seen as follows. Each ijjf has its support approximately concen- 
trated on configurations that display a particular macroscopic state of the pointer. This 
implies that the tp£ have negligible overlap and that the beable configuration (xj, . . . , x£) 
will be in the support of only one of the wave functions ipf, say As such it displays 
the kth macroscopic state of the pointer and hence the corresponding measurement out- 
come. Since in quantum equilibrium, the probability for the configuration (x", . . . , x£) 
to be in the support of ip^ is given by |cfc| 2 , also the quantum probabilities are recovered. 
Note that the wave functions ipf might be overlapping in which case the configuration 
of the system itself does not display the measurement outcome. However since measure- 
ment results are revealed to us be means of macroscopic records, this does not pose a 
problem. 

Further, if the different wave functions i/)*i/)£ stay approximately non-overlapping at 
later times, which is usually guaranteed by decoherence effects, then as one can easily 
verify from the guidance equations, only one of the wave functions, namely ipf.ipl, will 
play a significant role in determining the velocity field of the actual beable configuration 
of system and apparatus. Since the other wave functions ip^ipf, % 7^ k play only a 
negligible role, they can be ignored in the future description of the behaviour of the 
beable configuration. This is called an effective collapse. The effective collapse explains 
the success of the ordinary collapse in standard quantum theory. Unlike the collapse in 
standard quantum theory the effective collapse is not one of the axioms of the theory, 
but a consequence of the theory. 5 

Note that the property that wave functions of macroscopically different states are 
(approximately) non-overlapping in configuration space is a special property of the posi- 
tion representation. The same wave functions might not be non-overlapping in an other 
representation. This is very important for the choice of beable since it is unclear how 

5 In the de Broglie-Bohm theory, the wave function of a subsystem of the universe can be defined 
as the conditional wave function [4]. This wave function actually undergoes what could be called an 
actual collapse. 



16 



a pilot-wave model in which such wave functions are overlapping in the corresponding 
configuration space may reproduce the quantum predictions. Similar issues will arise 
when introducing field beables for quantum field theory, i.e. merely introducing some 
field beable which is distributed according to a quantum equilibrium distribution will 
not guarantee that the pilot-wave model reproduces the standard quantum predictions. 
In section 9.1.3 a non-trivial example of such a theory is considered. 

2.2 Field theories 

So far we have only considered systems that can be described by a finite number of 
degrees of freedom. The transition to a system that is described by a continuum of 
degrees of freedom is straightforward. The transition can be thought of as a replacement 
of the discrete label n of the coordinates q n by a continuum label x, i.e. q n — > 0(x). The 
sums over the discrete label n are then replaced by integrals over x and the derivatives 
d/dq n are replaced by functional derivatives 5/5<j)(x.). The corresponding field velocities 
d<p(t, x)/<9t will be denoted as 0(t, x). The field momenta 5L/5<j)(x.) will be denoted as 
n^(x). Usually we deal with a number of fields so that the fields also carry an extra 
discrete label r, i.e. r (x). 

The Hamiltonian formulation of a field theory involves some subtleties that are not 
present for a system that can be described by a finite number of degrees of freedom. 
However these subtleties are not important for the field theories we consider here. The 
interested reader is referred to [41-44]. 

2.2.1 Introducing field beables 

For the field theories we consider here, the constraints are of the form x m (0(x), IT^x)), 
m = 1, . . . , M. So for each m there is an infinite number of constraints, one for each 
point x in space (however, for convenience, we will often say that we have M constraints, 
instead of saying that we have an infinite number of them). The physical Hamiltonians 
will be of the form 

H = \H J ^ 3 yn^(x)/ i f s (x,y)n^( y ) + y(0), (31) 

where the superscript in hf s (x,y) denotes that this density can be a functional of 4>. 

On a formal level, the corresponding quantum field theory is obtained as follows. 
The canonical variables are replaced by operators and the commutation relations 

[0 r (x),fl 0s (y)] = i5 rs 5(x-y) (32) 

are imposed. These commutation relations are realized by the representation 

r (x) -> r (x) , IV (x) -> -i^—r^ • (33) 
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These operators act on wave functional *&((f>), for which the inner product is given by 

(*i|* 2 > = J {ll V A *l(<f>)M<f>) , (34) 

where T><fi r . = FJ X rf0 r (x). By introducing the field basis \cf>), where the \(f>) are eigenstates 
of the operators (f> r , i.e. r (x)|0) = r (x)|0), the wave functional \I/(0, t) can be regarded 
as the expansion coefficient (<fi\^(t)) of the state \^(t)). The dynamics of the wave 
functional is determined by the functional Schrodinger equation 

The representation for the field operators is the functional Schrodinger representa- 
tion. It is the natural generalization of the familiar representation (20) for ordinary 
quantum mechanics, but is not so widely used as for example the Fock representation. 
An introduction to the functional Schrodinger picture can be found in [56]. See also [57] 
for a detailed treatment of quantum field theory in the functional Schrodinger picture. 

A pilot-wave model can be constructed by introducing field beables r (x) whose 
dynamics is determined by the guidance equations 

M*) = J2[ d3 y h ^y)jJT-^ (36) 

,J °<Ps{y) 

where * = |\l/|e l5 was used. The quantum equilibrium distribution is given by |^(0, t) | 2 Y[ r T><p r . 



2.2.2 Regularization, renormalization and probability interpretation 

The presentation of the pilot-wave model was merely formal, lacking mathematical 
rigour. In particular, the measure rj r P0 r that appears in the inner product (34) and 
in the quantum equilibrium distribution is ill-defined. It was treated as an infinite- 
dimensional generalization of the Lebesgue measure, but such measures do not exist. 
Instead, one should adopt the infinite-dimensional analogue of a weighted measure [58- 
60]. With such a measure, the functional Schrodinger representation (33) needs to be 
modified in order to ensure that the operators are Hermitian. In addition, care is re- 
quired in the formulation of the Schrodinger equation, since such a measure is generally 
concentrated on distributions that are not smooth functions [61]. 

Another approach to make the Schrodinger picture well-defined was presented in [62- 
66]. There it was shown, at least for specific theories, that the functional Schrodinger 
picture can be made well-defined by a suitable regularization and renormalization. 

Given a well-defined functional Schrodinger equation, we also expect to be able to 
find well-defined guidance equations. 

The difficulties of course disappear when a regularization is introduced that makes 
the number of degrees of freedom finite. One such regularization consists in assuming 
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a finite spatial volume, together with some appropriate boundary conditions, and an 
ultra-violet momentum cut-off. The assumption of a finite spatial volume discretizes 
the possible momenta. The cut-off then further makes the total number of possible 
momenta finite. More explicitly, these assumptions imply that the fields r (x) have 
Fourier expansions 



where the possible momenta k are given by (rii, ri2, n^ir/L, with the n« G N such that 
|k| ^ A. L is the length of a side of the cubic spatial volume and A is the cut-off. In 
this way, the theory can be formulated in terms of the variables r k, of which there are 
finitely many. The Schrodinger equation will then be an ordinary wave equation for wave 
functions ^((/vki, . . . , r k n , t) and the beables will be a set of variables r k;, i — 1, . . . ,n, 
which are distributed according to (^(^rki, • • • ,0rk n ,^)| 2 in quantum equilibrium. The 
beables 4>rki will still give rise to a field </v(x) in physical space, through equation (37). 

This regularization makes both the Schrodinger equation and the guidance equation 
cut-off dependent. Removing the cut-off, by applying a suitable renormalization scheme 
and taking the cut-off to infinity is in general very hard (already on the level of the wave 
equation). On the other hand, it may be that nature provides a natural cut-off, so that 
the limit process is not required. 

Another possible regularization which makes the total number of degrees of freedom 
finite consists in assuming a finite spatial lattice. The degrees of freedom are then the 
"field" values at each lattice point. This type of regularization is perhaps more elegant 
than the cut-off regularization, since lattice gauge theories can be formulated in way 
which makes the gauge structure very similar to the gauge structure of the continuum 
theory (see for example the Kogut-Susskind approach [67, 68], in which space is treated 
as discrete and time as continuous). 

In case the regulators are not removed by some renormalization scheme, the pilot- 
wave ontology in the case of a lattice regularization would be different from the one in 
the case of the cut-off regularization that was described above (even though the lattice 
spacing actually also yields an ultra-violet cut-off, which is of the order of the inverse of 
the lattice spacing). In the case of a lattice regularization, the beables would represent 
degrees of freedom at the lattice points, whereas with the cut-off regularization the 
beables would be fields on space. 

In this paper, we maintain the formal approach, i.e. the approach outlined in the 
previous section, in developing the pilot-wave models. Although we have not studied in 
detail the effects of a more rigorous approach, we do not expect major deviations from 
the more formal formulation. 

2.2.3 Reproducing the quantum predictions 

As mentioned before in the context of the non-relativistic pilot-wave theory of de Broglie 
and Bohm, cf. section 2.1.3, assuming some beable that is distributed according to quan- 
tum equilibrium is not sufficient to guarantee that a pilot-wave model reproduces the 




(37) 
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predictions of standard quantum theory. The beables should yield, on the macroscopic 
level, an image of macroscopic objects. Therefore, it is important to consider wave func- 
tionals representing macroscopic states and find out what the typical field configurations 
are in quantum equilibrium. In particular, wave functionals representing different macro- 
scopic states should have their support approximately concentrated on different regions 
in field space. For example, whether Schrodinger's cat is alive or dead, or whether there 
is a cat in the first place, depends on the values the field beables take and if the wave 
functionals representing a live cat and a dead cat would be very much overlapping, then 
there would be nothing in the field beable revealing whether we are dealing with a live 
cat or a dead cat. Similarly, in measurement situations, the beables should display a 
definite state of the macroscopic pointer. If the wave functionals corresponding to dif- 
ferent states of the macroscopic pointer have significant overlap then there is no hope 
that the pilot-wave model wil reproduce the standard quantum predictions. 

In [69] Saunders expressed a worry for using fields as beables. He expressed some 
doubts whether localized macroscopic bodies are represented by localized field beables. 
If not, he claimed, it would be unclear how a pilot-wave model with field beables repro- 
duces the quantum predictions. However, although it is true that a pilot-wave model 
in which localized macroscopic bodies are represented by localized fields may reproduce 
the quantum predictions, it is by no means a necessary requirement. Macroscopic bodies 
could also be displayed by non- localized fields. For example macroscopic objects could 
be represented by fields that differ in magnitude in certain region of physical space. But 
although a field configuration representing a localized macroscopic object need not be 
localized, it is still desired that the field configuration displays the image of localized 
macroscopic object locally, i.e. it is desired that the presence of the macroscopic object 
can be inferred from the characteristics of the field configuration in the spatial region 
where the object is supposed to be located. Although an ontology which does not satisfy 
this property can possibly be maintained in some cases, it would be rather far removed 
from our everyday experience of the macroscopic world. 

2.2.4 Non-locality and Lorentz invariance 

We want to end this section with a note on non-locality and Lorentz invariance. The 
pilot-wave models we present here are non-local. Because of Bell's theorem the non- 
locality is unavoidable. The pilot-wave models are further formulated with respect to a 
preferred frame of reference. In the context of a relativistic quantum field theory this 
will imply pilot- wave models that are not Lorentz invariant. Nevertheless, in quantum 
quilibrium these pilot-wave models will reproduce the standard quantum theoretical 
predictions. Therefore the empirical predictions of the pilot-wave models will be Lorentz 
invariant, at least when the quantum theoretical predictions are Lorentz invariant [16, 
17,70,71]. There are some attempts to formulate Lorentz covariant pilot-wave models. 
An overview of these is given by Tumulka [38]. See however also Valentini [72] for an 
argument that the fundamental symmetry of pilot-wave theories is Aristotelian rather 
than Lorentzian. 
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3 Schrodinger field 



The first field theory that is considered is the quantized non-relativistic free Schrodinger 
field. The Schrodinger field can be quantized by imposing either bosonic or fermionic 
commutation rules. Here we only discuss bosonic quantization. Fermionic quantum field 
theories will be discussed in section 9. 

The pilot-wave approach to the bosonically quantized Schrodinger field was discussed 
before by Takabayasi [21] and Holland [22] and [17, pp. 449-451]. Our presentation 
differs from that of Takabayasi and Holland. We also obtain a slightly different pilot- 
wave model. 



3.1 Pilot-wave model 

First consider the quantization of the Schrodinger field. The classical field equation is 
just the non-relativistic Schrodinger equation 

which can be derived from the Lagrangian 

L = Jd 3 x (>> - ^*) + ^* V 2 ^) . (39) 

In order to pass to the Hamiltonian formulation it is convenient to express in 
terms of real and imaginary parts. 6 We write ip — ip r + r0j. The canonically conjugate 
momenta, the Hamiltonian and the constraints are respectively given by 

SL S L 

IV = — j- = ipi , = —r- = -ip r , (40) 

H = J d 3 x fyvVVr + A^ 2 A) , (41) 

xi = -ipi, X2 = + ip r . (42) 

Since [xi( x )> X2(y)]p = — 25(x — y), the constraints are second class. 

As explained in section 2.1.1, we proceed by separating the unconstrained variables 
from the constraints by performing a canonical transformation. Such a canonical trans- 
formation is given by 

Vv = ^(0 + 0O , n^ = — (n^ + n^) , 

A = 72 (n ^ " n0,) ' n ^ = 72 (0 ' " 0) ' (43) 



6 See [73] for an overview of possible derivations of the Hamiltonian formulation for the Schrodinger 
field. 
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In terms of the new canonical variables 0, 11^, 0' and Up, the constraints read 0' = U^r = 
0, so that and 11^ are unconstrained variables. The Hamiltonian for the unconstrained 
variables reads 

H = -^-Jd 3 x (n^V 2 n + 0V 2 0) • (44) 
The corresponding Hamilton equations reduce to the second order differential equation 

0+^0 = 0. (45) 

Note that this equation of motion corresponds to the square of the non-relativistic 
dispersion relation. 

Quantization proceeds by associating field operators and II^ to the canonical vari- 
ables and by imposing the standard commutation relations 

[0(x),n^(y)] = W(x-y). (46) 
The corresponding Hamiltonian operator is given by 

H = — — — f d?x (n V% + 0V 2 0) . (47) 



4m _ 

Using the functional Schrodinger representation 



0^0, n ^-i^, (48) 

the functional Schrodinger equation 

if = / (f x j* y (-i^Mx,,)^ + ^(x)Mx,y)0(y)) * (49) 
is obtained, where 

/ i (x,y) = --^V 2 5(x-y). (50) 
2m 

A field beable can now be introduced, which evolves according to 

^).| AA(x>y) _^__L v2 ^|_ (51) 

and for which the quantum equilibrium density is given by |\l/(0, t)| 2 . 

3.2 Comparison with alternative formulations 
3.2.1 Alternative Hamiltonian formulation 

Usually the quantized Schrodinger field is formulated in terms of the field operators ip 
and ■0*, which satisfy the commutation relations 

[#(x),^(y)] = tf(x-y). (52) 



22 



A systematic way of deriving these commutation relations involves constructing the 
Dirac bracket for the classical fields -0 and -0* (which are treated as independent fields), 
and their corresponding momenta 11^ and H^* [73], which are given by 

n^ = r0*/2, IV = -r0/2. (53) 

The formulation in terms of the operators -0 and -0* is equivalent to the formulation 
in terms of and 11^. The equivalence is obtained by using the relations -0 = -0 r + vpi 
and -0* = -0 r — iipi, by using the operator form of the canonical transformation (43), 
i.e. ip r — (0 + 00/^/2 and ipi = (IT^ — T\.@)/y/2, and by imposing the constraints as 
operator identities, i.e. 0' = Up = 0. One can easily see that these relations imply the 
equivalence of the commutation relations (46) and (52). 

In terms of the field operators -0 and -0* the Hamiltonian operator (47) has the 
familiar form 

H = — — / d 3 x$*V 2 $ + E , (54) 
2m J 

where 

E = ^-jd 3 xj d 3 yV x 5(x - y) • V x 5(x - y) (55) 

is an infinite vacuum energy. It is common practice to drop the infinite constant. This is 
usually done by assuming a suitable operator ordering. Here we will leave the operator 
ordering as it is. The removal of the infinite constant in the Hamiltonian operator for ip 
and -0* would imply the addition of this infinite constant in the Hamiltonian operator 
for (f) and n^. The appearance of this divergent term clearly shows the need for a 
regularization. 

So our formulation of the quantized Schrodinger field is equivalent to the more fa- 
miliar formulation in terms of the field operators -0 and -0*. In retrospect, we could 
have developed the pilot-wave model starting from this formulation, by realizing the 
commutation relations of the operators -0 and -0* by the representation 

(*"£)■ (56) 

3.2.2 Alternative pilot-wave approaches 

Note that by introducing a field beable 0(x) for the field operator 0(x) we have in fact 
introduced a beable for the operator 0v(x), namely the field ip r = (p/y/2. Alternatively, 
we could have introduced a field beable n^(x) for the momentum field operator n^x), 
which is related to the operator -0j(x). This underdetermination can be traced back 
to the underdetermination in the canonical transformation that separates the uncon- 
strained degrees of freedom from the constraints, cf. section 2.1.1, since one can always 
perform a canonical transformation that interchanges the role of fields and momenta. 
Similar underdeterminations were pointed out for the case of the electromagnetic field 
[74]. 
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There is also a way to introduce beables for fa and fa simultaneously (instead of for 
either one of the two), and hence for if>. For example, starting from the pilot- wave model 
presented above, with beable fa = <$/ a/2 corresponding to fa, one could also introduce a 
beable fa corresponding to fa, by means of Holland's local expectation value, cf. section 
2.1.2. In the field representation (56) the operator fa is represented by the operator 
-j^jqi so that the corresponding beable is given by 



fa{x.,t) = Re 



TO,*)! 1 



1 6S(<f),t) 



4>=<t>(x,t) 



y/2 50(x) 



, (57) 

0=0(x,t) 



where the expressions on the right hand side are evaluated for the actual value of the 
field beable 0(x, t). The beable ip for the operator ip is then given by 

1 { ss\ 

$ = fa + ifa = — \6 + i—j , (58) 

where the right hand side is evaluated for the actual beable configuration <fi. 

Another way to obtain a beable if) for if) would be to first introduce a beable 11^ for 
the momentum field operator 11^ in the standard way (using the continuity equation 
to obtain a guidance equation), and hence for fa, and to subsequently introduce a 
additional beable for fa by means of Holland's local expectation value. The pilot-wave 
model as such obtained would be different from any of the models discussed above. 

In the approach of Takabayasi [21] and Holland [22] and [17, pp. 449-451], the beable 
is also given by a complex field if)(x), with a form similar to that of (58). The difference 
merely arises from the fact that Takabayasi and Holland use a representation of the field 
operators that is different from (56). 

All of the possible pilot-wave approaches just described concern a field in physical 
space, so that they provide a clear ontology and hence there is no a priori reason to prefer 
one ontology above the other. Of course, it could be that either ontologies fails to yield 
an image of macroscopic objects and hence fails to reproduce the standard quantum 
predictions and could be excluded on this basis. Note also that theories with different 
ontologies that reproduce the standard quantum predictions are evidently empirically 
indistinguishable. In quantum non-equilibrium these theories do not yield the standard 
quantum predictions and they might under these circumstances be distinguishable. 

However, the Schrodinger field describes identical bosonic particles, 7 so that, despite 
the many applications of the theory, it can not really faithfully describe macroscopic 
objects, like pointer needles; since such objects are composed out of fermions, they 
should be described by a fermionic field theory. Therefore we can not really discuss 
how our pilot-wave model in terms of the field beable is able to reproduce the quantum 
predictions. We return to this in section 3.4. 



7 Note that in writing about particles here, we have employed the ordinary parlance of standard 
quantum theory; we are not referring to actual point-particles. Throughout the paper, it should always 
be clear from the context what meaning we have in mind. 
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3.3 Valentini's suggestion for dealing with constraints 



3.3.1 Dirac field 

In [23,24] Valentini tried to develop a pilot-wave model for the quantized Dirac field 
(which describes relativistic spin-1/2 particles). With the quantization of the Dirac 
field constraints appear that are similar to the ones we encountered in the case of the 
Schrodinger field. It is interesting to consider how Valentini dealt with the obstacle 
of constraints, because it is very similar in spirit to the approach adopted here. In 
fact, since Valentini was dealing with fermionic quantization, he was facing additional 
obstacles in the construction of a pilot-wave model. However, these additional obstacles 
will be discussed in detail in section 9.2. 
For the Dirac field the constraints read 

H^ = i^/2 , H c = -i^/2 . (59) 

Note that these constraints are completely similar to the constraints (53) for the Schro- 
dinger field; there is just an extra spinor index in this case. Valentini realized that it 
would make no sense to impose the guidance equations 

-T = ^' (60) 

since rather than telling us the rate of change of the field beable, the functional derivative 
5S/5ip a would tell us the actual value of the field beable. Valentini reasoned that the 
appearance of the constraints, and hence the reason for the difficulties, was due to the 
fact that the Dirac equation is first order in time derivatives. As a solution, Valentini 
proposed to use a second order formulation of the theory instead. The second order 
equivalent of the first order Dirac theory is the Van der Waerden theory. Since the 
Hamiltonian formulation of the Van der Waerden theory does not involve constraints, 
the above problem disappears. 

Let us comment on this. First of all, we want to say a word on the particular choice 
of the guidance equation (60). Valentini obtained the guidance equation by equating 
5S/5<f) r and H r {4>, <fi), where S is the phase of the wave functional and where H r ((f>,<f>) 
are the classical momenta conjugate to the fields <f> r . This prescription is inspired by the 
Hamilton- Jacobi formulation of Hamiltonian systems. In a number of familiar situations, 
like e.g. in the non-relativistic pilot-wave theory of de Broglie and Bohm, this prescrip- 
tion yields the same guidance equation as the one that can be derived by considering 
the continuity equation for the probability density |\I/| 2 . However, this prescription does 
not work in general. The Dirac field is just one example where it fails. In Appendix 
A, another simple example is given which, unlike the Dirac theory, involves bosonic 
quantization and has a Hamiltonian formulation that does not involve constraints. 

Despite of this, Valentini's insight of using a second order formulation of the theory 
is in agreement with our approach. In our approach the constraints were separated from 
the unconstrained variables and only the latter were quantized. In effect, the first order 
Schrodinger equation was thereby replaced by a second order equation (45) for the field 
by 0, which is in accordance with Valentini's suggestion. 
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3.3.2 String field theory 



In [75] Weingard tried to develop a pilot-wave model for the second quantized bosonic 
string (see also [76] for a summary). In doing so, Weingard also faced the problem of 
constraints. Weingard adopted the line of reasoning of Valentini and suggested to look 
for a second order reformulation of the theory. Such a second order reformulation was 
unknown to Weingard, which forced him to halt the attempt to develop a pilot-wave 
model. However, the theory of the second quantized bosonic string is formally similar 
to the quantized Schrodinger field and therefore a pilot-wave model can be developed in 
a similar way as for the Schrodinger field. Let us sketch how this would go. 

First recall the basic ingredients of string theory in the light-cone gauge. A clas- 
sical string in (D — l)-dimensional space traces out a two-dimensional surface in D- 
dimensional space-time, called the world-sheet. The world-sheet can be represented by 
the coordinate functions x^(a, r). It can also be represented in terms of light-cone co- 
ordinates with i = 1, . . . , D — 2, where x + and x are coordinates along the 
light-cone and x % the D — 2 remaining spatial coordinates. In the light-cone gauge the 
parameter r is proportional x + and the coordinate function x~ can be expressed in terms 
of the x l and a constant of the motion Xq . In the resulting Hamiltonian formulation, 
the parameter r plays the role of time and the dynamical degrees of freedom are the 
functions x l {a). 

In the first quantized theory, the states are functionals ip(x l (a) : r) [57] and a pilot- 
wave model can be constructed by introducing beables x l (a) which represent a curve in 
D — 2 dimensions [75]. (Hereby the degree of freedom Xq has been ignored. We leave 
aside the question whether or not an additional beable should be introduced for x^ or 

x~ .) 

In the second quantized theory, one has operators ip(x l (cr)) and i/;*(x l (cr)), which 
satisfy the commutation relations 

Mx),r(y)] = 6(x-y) (61) 

and the Hamiltonian is of the form 

H = J VxVy$*(x)h(x,y)$(y) , (62) 

where h(x,y) is some kernel [75]. These equations are formally similar to the ones for 
the quantized Schrodinger field, cf. (52) and (54), and hence a pilot-wave model can be 
developed in a similar way as for the Schrodinger field. Without going into detail, this 
would result in a beable 0(x*,r) that is being guided by a wave functional \I/(0(:r l ), r). 
The beable <f)(x l , r) would then be a real-valued functional, defined on the space of curves 
x l (a) in a (D — 2)-dimensional space. Whether such a beable is capable of providing 
an image of macroscopic objects (so that, in particular, measurement results would be 
recorded and displayed in the beable) is of course an open question. 
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3.4 Some elementary states 

As explained in section 2.2.3, in order to find out whether a pilot-wave model is capa- 
ble of reproducing the standard quantum predictions, it is important to consider which 
fields make up the support of wave functionals representing macroscopic objects. In 
particular, wave functionals representing macroscopically distinct states should have 
approximately disjoint supports. However, as mentioned before, the quantized Schro- 
dinger field describes identical bosonic particles and can therefore not really faithfully 
describe macroscopic objects like pointer needles. Hence we can not really discuss how 
the pilot-wave model is able to reproduce the quantum predictions. Nevertheless, it is 
still interesting to look at some familiar states and consider the asociated field distribu- 
tions. This is because similar results apply to other bosonic field theories discussed in 
this paper, like for example the electromagnetic field. 

We will subsequently consider the vacuum state, coherent states and states which 
describe a system with a definite number of particles. Since we are merely interested 
in the associated field distributions we will not give the associated guidance equations. 
Some discussion of the guidance equations corresponding to these states, in the context 
of the scalar field and the electromagnetic field, can be found in respectively [15-20] and 
[10]. 

3.4.1 Vacuum state 

As is well known the operators '0*( x ) an d V'(x) respectively create and annihilate a 
particle at the position x. The operators a^(k) and a(k), which are introduced through 
the relations 

£(x) = ^ f dV-a(k) , </>(x) = ^ / ^e- ik *at(k) , (63) 

respectively create and annihilate a particle with momentum k. The vacuum, which 
contains no particles, is represented by the state |\E'o(^))- It satisfies V ; (x)|\[ / o(^)) — for 
all x or equivalently a(k)\^ (t)) = for all k. 

In the functional Schrodinger picture, we find by using (56) that the vacuum wave 
functional *o(0, t) = (<f)\^o(t)} is determined by 

(* (x) + ih )*° (M = - (64) 

The solution reads 

MM = iVexp (-± J rf 3 x0(x) 2 - iE t) , (65) 

where N = rix 71 " -1 ^ 4 i s an infinite normalization constant. It can easily be verified that 
ty ((j),t) satisfies the functional Schrodinger equation. The ground state is a Gaussian 
wave functional which is centered around the field configuration = 0. 
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3.4.2 Coherent state 



Another state of interest is the coherent state \^ a (t)) [77]. The coherent state \ty a (t)) 
satisfies 

a(k)|* a (t)) =a(k)e- ik2t / 2m \^ a (t)). (66) 
In case we can take the Fourier transform of a(k), we also have 

$(x)\* a (t))=a(x,t)\* a (t)), (67) 

where 

a(x,t) = y rf 3 A;e ik - x - ifc2 */ 2m a(k) . (68) 

In the functional Schrodinger picture, the coherent state ty a ((f),t) = (<p\^/ a (t)) satis- 
fies 

71 ( 0(x) + s^j) t} = a(x ' t} * a(0, t} ' (69) 

The solution is given by 

* a (0,t) = iVexp ^-i y rf 3 x (V(x) - ^(x.t)) 2 - i£b^ , (70) 

where N = (Ylx 71 ^ 1 ^) ex p( — / d 3 xaf) is a normalization constant and where a = a r + 
iojj was used, with a r and aj real. 
Writing 

*a(<M) = (^11^ V4 ^) ex P ( - ^ / ^ (^W - v^a r (x,0) 2 

+ iv^2 y rf 3 x (<£(x) - v^a^x, «i(x, t) - iE ?j , (71) 

we see that the field density |^/ Q (0, t) | 2 is a Gaussian centered around the field A/2a r (x, t). 
In the case the function a is zero we have the ground state wave functional. 

It is clear that it is sufficient that the real parts of cci and a 2 differ significantly 
in just a region of physical space in order for the wave functionals \I/ ai and ty a2 to 
be non-overlapping. Nevertheless, these wave functionals do not describe macroscopic 
systems in non-relativistic quantum theory. Of course we are considering a theory which 
describes only identical bosonic particles, but there is an additional reason. The reason 
is that we think of a macroscopic system, like a pointer needle, as a system composed 
out of an approximately definite number of particles. To the contrary, the coherent state 
is a superposition of states with all possible particle numbers. On the other hand, if we 
consider electromagnetism, then coherent states for which the mean photon number is 
sufficiently large will describe familiar macroscopic states of light. Therefore the above 
discussion becomes important for the pilot-wave theory for the electromagnetic field 
(see section 4). States for the Schrodinger field describing a definite particle number are 
discussed in the following section. 
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3.4.3 State describing a definite number of particles 

An arbitrary state describing n particles is given by 

|^(t)) = -i=y d 3 Xl ...J rf 3 ^(x 1 ,...,x n ,t)^(x 1 )...^(x n )|^ (t)), (72) 
where the completely symmetric expansion coefficients ip{x.i, ■ ■ ■ , x n , t) are given by 

^( Xl , ...,x n ,t) = 4=(*o(OI^K) • • • ?(xx)|^(i)) . (73) 



The requirement that \^^{t)) satisfies the field theoretical Schrodinger equation implies 
that ^( x i! • • • ) x n 5 1) satisfies the non-relativistic Schrodinger equation. Consider the 
wave function ip{x.i, . . . , x n , t) to be normalized to one, so that the state \^^(t)) is 
normalized to one. 

For a single particle the wave functional reads 

*V>(<M) = j d 3 x^,t)^= (V(x) - V (<f>,t) 

= V2 J d 3 a^(x,t)0(x)^ o (<M) • (74) 

The wave functional corresponding to any number of particles can be written in a similar 
fashion. 

In the pilot-wave theory of de Broglie and Bohm for non-relativistic quantum theory, 
two states representing a single particle localized at different regions in physical space 
are non-overlapping. With the field ontology this does not seem true anymore. There 
is a strong indication that wave functionals describing a single particle will in general 
have significant overlap, regardless of whether these states describe localized particles. 8 
To see this consider the following. 

The field density corresponding to a one-particle state ^ is given by 

|^(0)| 2 = 2|a| 2 |vl/ o (0)| 2 , (75) 

where 

a = / rf 3 x^(x)0(x) . (76) 



The density reaches a maximum for fields that satisfy 5\^^\ 2 /5(f) = 0. This implies that 
either 9 

a = (77) 



8 Valentini, private communication. 

9 A similar calculation was done in [23] for the scalar field. 
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For a field <p f° r which a — 0, the density |\1/^(0)| 2 is zero and hence for such a field 
a minimum is reached. Maxima are reached for fields which satisfy (78), i.e. for fields 
which are a linear combination of ip and vp*. 

Consider now the superposition vpjv^i+ifo) = ^ (^Vi + ^^2)5 °f two one-particle 
states and \l/^ 2 , with iV a normalization factor. This state has maxima for fields 
that are linear combination of ipi + ip2 and ip\ + On the other hand, if we assume 
that and \l/^ 2 have negligible overlap, then 

l^ 1+ ^ 2 )r~iV 2 (l^il 2 + l^ 2 | 2 ) • (79) 

The right hand side reaches maxima for fields which are a linear combination of ipi and 
ip* and for fields which are a linear combination of t/j 2 and These do not correspond 
to the maxima obtained for \^/ N ^ 1+ ^\ 2 . Because of (79), this suggests that in general 
the states and will have significant overlap. 

For clarity, in the special case the wave functionals and \P^, 2 represent a particle 
localized in different regions in physical space, tpi and ip2 will have negligible overlap in 
physical space. However, as the above analysis shows, this does not imply that the wave 
functionals and have negligible overlap in field space. 

We think of a macroscopic system, like a measurement needle, not only as an approx- 
imately localized system, but also as a system composed of a large, but approximately 
fixed number of particles. Therefore it would actually be more interesting to consider 
wave functionals for such systems and see in which cases they might be non-overlapping. 
Recently some numerical simulations have been performed by Schmelzer which indicate 
a decrease in overlap that goes exponentially with the number of particles [78]. 

4 Free electromagnetic field 

We now turn to the free electromagnetic field. In this and the following sections, the 
components of four-vectors are labeled by Greek indices, e.g. /x and u, and the compo- 
nents of spatial three- vectors are labeled by Latin indices, e.g. % and j. For these indices 
we use Einstein's summation convention. 

We consider two pilot-wave approaches to the electromagnetic field. One is by Bohm 
and takes a tranverse potential as beable. The other one is by Valentini and takes an 
equivalence class of potentials as beable. While these ontologies are different, these 
models appear to be empirically equivalent, not only at the level of quantum equilibrium, 
but also for non-equilibrium distributions. 

4.1 Hamiltonian formulation 

The Lagrangian and the Euler-Lagrange equations of motion for the free classical elec- 
tromagnetic field, described by the vector potential = (A , A*), are respectively given 
by 

L = ~ I d 3 xF^F^ , = d»A v - d v A» , (80) 
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d^ u = 0. (81) 

The theory has a local symmetry 

N 1 -> A» - d^9 . (82) 

The arbitrariness of 9 as a function of space-time signals gauge freedom. 10 This gives 
rise to first class constraints in the Hamiltonian picture. 
Indeed, in the Hamiltonian picture we have [41-44]: 

SL 5L 

n ^ = 7T- = ' n Ai = 71 - = A + 9iA), (83) 

H = Jd 3 x Qn Ai n Ai + ^FijFij + AodiUA^j , Fij = diAj - djA, , 

Xi = H-a , X2 = diU Ai , (85) 

where the constraints Xi an d X2 are first class. Since II Ai = A\ + diA = —E^ with Ei 
the electric field, the second constraint can be written as diEi = 0, which is just one of 
the free Maxwell equations, namely the Gauss law. 



4.2 Bohm's approach: transverse potential as beable 

In presenting his pilot-wave model for the electromagnetic field, Bohm started with 
imposing the Coulomb gauge d{Ai = [3]. (From the classical equations of motion 
it then follows that V 2 A = 0, so that, assuming that A vanishes sufficiently fast at 
spatial infinity, A — 0.) Bohm then expressed the transverse parts of the fields in terms 
of Fourier modes, which form unconstrained canonical pairs. The pilot-wave model was 
then introduced for the quantum field theory that results from quantizing those Fourier 
modes. Belinfante [80, pp. 198-209] and Kaloyerou [9-11], who discuss Bohm's model 
in detail, adopted a similar approach. 

The quantum field theory for the transverse modes of the field can actually be ob- 
tained more directly by identifying unconstrained canonical pairs as gauge independent 
degrees of freedom and by quantizing these (as discussed in section 2. 1.1). 11 We follow 
this approach to present Bohm's model because it immediately makes clear that the 
model is gauge independent, not only at the statistical level in quantum equilibrium, 
but also at the fundamental level. This also immediately answers the questions raised 
by Kaloyerou [10, 11] concerning the gauge invariance of the pilot-wave model. 

10 It is straightforward to show that the transformations (82) are all gauge transformations in the 
sense defined in section 2.1.1. More details can be found in [79]. 

11 In view of the last paragraph of section 2.1.1, it is not surprising that the same theory can be 
obtained without gauge fixing. 



31 



4.2.1 Quantization of gauge independent degrees of freedom 

In order to quantize only the gauge independent degrees of freedom, the first step is to 
perform a canonical transformation so that in terms of the new canonical variables, the 
constraints are given by some of the momenta. Since the momentum Ua is constrained 
to be zero, A can immediately be identified as a gauge variable. In order to identify the 
remaining gauge variable, we write = Af + Af, where Af and Af are respectively 
the transverse and longitudinal part of A iy defined by 

Aj=^-^y j , Af = ^A,, (86) 
and where the integral operator V -2 is given by 

^=-h^\- < 87 » 

The fields Aj and IT^. have vanishing Poisson brackets with the constraints and hence 
they are gauge independent degrees of freedom (as is well known [41-44,81]). The field 
Af does not have vanishing Poisson brackets with the constraint X2 and is the remaining 
gauge degree of freedom. 

In order to write the constraint \2 = c^II^ and the gauge degree of freedom Af in 
terms of a canonical pair, the following canonical transformation of the variables A t and 
Tl-Ai is performed: 

n*(x) = (2^372 E / ^- ik ' x 4(k)n ?( (k), 

where the new (complex) variables <ft(k) and II 9i (k) form canonical pairs. We also 
introduced the (real) orthonormal polarization vectors e'(k), I — 1, 2, 3, which we choose 
to obey £ 3 (k) = hi/k and e'(k) = e l (— k) for I — 1, 2. From the last relation and the fact 
that Ai and U A . are real, it follows that qi(k) = g z *(— k) and n gi (k) = k) = II* (—k) 

for 1 — 1,2, with similar relations for q 3 and U q:j . 
From the relations 

A l = E / d^e\(k) qi (k) , Af = ^ / ^e^(k)g 3 (k) , 

2 

n i = (2^5/2 E / ^- ikx 4(k)iWk) , = | rf 3 A;e- ik -% 3 (k)n, 3 (k) , 

(89) 



we see that the transverse parts of the fields are described by the variables q±, qi, n gi , IIq 2 . 
The longitudinal parts of the fields are given by g 3 and II 93 . Because the constraint 
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dillAi = implies 11^ = 0, it reduces to U q3 = in terms of the new variables, so that 
g3 is the remaining gauge variable. 

The Hamiltonian for the gauge independent variables qi,q2, n gi and U q2 is given by 



2 



(90) 



As is well known, this resulting theory is readily obtained by imposing the Coulomb 
gauge and by expressing the transverse fields in terms of Fourier modes. 
After quantization, we obtain the functional Schrodinger equation 

for the wavefunctional ty(qi,q2,t). 
4.2.2 Pilot-wave model 

Beables can now be introduced in the way outlined in section 2.2.1. They form a pair 
of fields (51,52) that evolve according to the guidance equations 

« k) ^- <^ 2 - (92) 

The quantum equilibrium density for the field beables is given by \^(qi, <?2, t)\ 2 . 

The associated ontology in space is given by a transverse vector field Aj (x) that 
is obtained from the fields (gi(k), g 2 (k)) using the relations (89). The dynamics can 
be expressed directly in terms of the transverse field Aj (x) as follows. The functional 
derivative 

can be decomposed into a transverse and a longitudinal part 

M,(x)- (2^)3/^7 £lW ^(k)' 

1 / 'd 3 A;e- ik - x e?(k)— 4^- (94) 
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^f(x) 




V 2 
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<L4f(x) 


V 2 





(27T)3/2y ^ ^g 3 ( k ) 

By inverting the first relation the functional derivatives 5/5qi(\t), I = 1,2, can be written 
in terms of S/5AJ (x). Expressing also the fields qi(k), I — 1, 2, in terms of the transverse 
vector field (x), by inverting the relations (89), the functional Schrodinger equation 
(91) can be rewritten as an equation for the wave functional ^(Af,t) 
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The guidance equation (92) then reads 




5S 



(96) 



5AJ 




and the quantum equilibrium distribution \^(A[ ,t)\ 2 VA T (where VA T ~ Vq{Dq 2 ). 

The ontology in Bohm's pilot-wave model is given by the field Aj. Just as in the 
case of the Schrodinger field, cf. section 3.2.2, and as noted by Baumann [74], alternative 
choices are possible. One could for example consider an ontology that is closer to that of 
classical electromagnetism in terms of the electric and magnetic field. For example one 
could use Holland's local expectation values to introduce a beable B i = e^kdjA^ for the 
magnetic field operator Bi and a beable Ei = —5S/5AJ for the electric field operator 



4.2.3 Macroscopic states 

Wave functionals representing macroscopically distinct states of the electromagnetic 
field describe distinct classical electric and magnetic fields. This implies in particular 
that the wave functionals give approximately disjoint magnetic field distributions. Since 
different magnetic fields Bi correspond to different fields Aj, we have that these wave 
functionals also give approximately disjoint distributions for the fields Aj . 

This feature is important since it implies that, in the context of a more complete 
theory which also includes a description of matter, like for example quantum electro- 
dynamics, measurement results can get recorded and displayed in the field beable Af. 
This observation was used in [26,27] to argue that it is actually sufficient to have only 
beables for the electromagnetic field in order to have a pilot- wave model that reproduces 
the predictions of, say, quantum electrodynamics. We discuss this approach further in 
section 10. 

In section 3.4, we considered some particular wave functionals for the Schrodinger 
field and the corresponding field distributions. The results of that discussion also apply 
to the case of the electromagnetic field. Wave functionals corresponding to a single 
photon will in general have significant overlap. On the other hand, coherent states that 
are "sufficiently distinct" are non-overlapping. With "sufficiently distinct" we mean for 
example that the states correspond to a different average photon number or that they 
correspond to photons with different momenta. For wave functionals corresponding to a 
definite but high number of photons it is unclear whether they can be non-overlapping. 
Note that a macroscopic state for the electromagnetic field is not a state with a definite 
number of particles (whereas we consider macroscopic states for ordinary matter to be 
states with an approximately definite number of particles), since for these states the 
expectation values of electric and magnetic field operators are zero. Instead a macro- 
scopic state contains, just as coherent states (apart from the ground state), an indefinite 
number of particles. 
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4.3 Valentini's approach: equivalence class of potentials as be- 
able 

Valentini has a different approach to gauge theories [23-25]. But while the ontology is 
different than that in Bohm's theory, the theories appear to be empirically equivalent, 
even in non-equilibrium. 

Valentini's pilot-wave model is best understood in the context of the quantization 
procedure in which constraints are imposed as conditions on states. Therefore we start 
with recalling this quantization procedure for the electromagnetic field. 



4.3.1 Quantization by imposing the constraints as conditions on states 

Instead of quantizing only gauge independent degrees of freedom, the electromagnetic 
field can also be quantized as follows (cf. section 2.1.1). All canonical variables are 
quantized as if there were no constraints, i.e. 

[A (x),n Ao (y)] = W(x - y) , [A(x),fU.(y)] = i<M( x " y) , (97) 

where the other fundamental commutation relations are zero, and the constraints are 
imposed as conditions on states, i.e. 

£i|*) = fU|tt) = 0, X2\*) = d i nA i \*) = 0. (98) 

Because the operators satisfy the standard commutation relations, the functional Schro- 
dinger representation can be used: 

a (x)^a (x), n Ao (x)^-i 



<L4 (x) ' 

A(x)^A(x), n Ai (x)^-i^-^. (99) 

The corresponding functional Schrodinger equation for the wavefunctional ty(Ao, Ai, t) = 
(A ,Ai\ty(t)} reads 

r J3 / i 5 2 i 



^-] dx {--2jAM^A F ^r (100) 

and physical states further have to satisfy the constraints (98), i.e. 

u-r ' (101) 

Sty 

a,- = o. (102) 

The quantum field theory so obtained is equivalent to the quantum field theory that 
was obtained by quantizing only the gauge independent degrees of freedom. This can 
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be seen as follows. The first constraint (101) simply states that \& does not depend on 
A . The second constraint (102) further implies that 

n mS * 0, (103) 



SAf V 2 6Aj 

i.e. physical states should also not depend on the longitudinal part of the vector potential. 
Hence physical states are of the form ty(Af,t). For these states the functional Schro- 
dinger equation (100) reduces to 

■£ 4M-*Ar w h (104) 

which is the Schrodinger equation (95) that was obtained by quantizing only the gauge 
independent degrees of freedom. 

In order to obtain full equivalence of the quantum field theories, one also needs to 
show the equivalence of the inner product. As mentioned in section 2.1.1, the definition 
of the inner product requires some care in the case constraints are imposed as conditions 
on states. In the Hilbert space of unconstrained wave functionals ^/(A ,Ai), one can 
define the inner product as 

(*ii* 2 ) = j va (n™^ ( ios ) 

However, physical states do not depend on A and A\ and therefore for physical states 
the expression (^1^2) becomes proportional to the infinite factor J T>A T>A L , the so- 
called "gauge volume" , which makes the wave functionals non-normalizable. The solu- 
tion is to factor out this infinite gauge volume from the expression (^1^2) and to define 
the finite remnant as the inner product for physical states, i.e. 

(*i|*2> P ft = J VA T *l* 2 . (106) 

This inner product is the same as the one used in the case where only the gauge inde- 
pendent degrees of freedom of the electromagnetic field are quantized. 

In general the gauge volume can be factored out by using the Faddeev-Popov proce- 
dure. This procedure basically consists in inserting the quantity 5(f(A , Ai))A F p,f(A , A*) 
in the integral in (105), where the delta function imposes a gauge fixing (i.e. the condi- 
tion f(Ao, Ai) = selects a unique representant from each set of gauge equivalent fields) 
and where Appj is the corresponding Faddeev-Popov determinant which ensures that 
this procedure is independent of the choice of gauge. In the case of the electromagnetic 
field there are many suitable gauge fixings. For example the form (106) for the inner 
product is readily obtained by choosing the gauge d{Ai = 0, A = 0. In the case of 
non-Abelian Yang-Mills theories it is more difficult to apply the Faddeev-Popov pro- 
cedure, because a large class of possible gauge fixings is ruled out [82-84]. Therefore, 
for non-Abelian Yang-Mills theories, the Faddeev-Popov procedure is often applied only 
locally, which is fine within a perturbative treatment. 
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4.3.2 Pilot-wave model 



Valentini did not really follow the above quantization scheme to arrive at his pilot-wave 
model. Instead he started by arguing that the degree of freedom Aq in the theory of the 
classical electromagnetic field is a mathematical artifact, arising from the insistence on 
Lorentz invariance. Valentini therefore dropped the degree of freedom A from the out- 
set. In the resulting quantum theory, the quantum state \I/ is a wave functional defined 
on the space of vector fields Ai which satisfies the functional Schrodinger equation 

'W A HsSs+i**)'- (107) 

together with the constraint 

9,^ = 0. (108) 

Valentini obtained the constraint by requiring that the wave functional be invariant 
under infinitesimal time-independent gauge transformations of the form 

Ai-tAi + diO. (109) 

This quantum field theory can also be obtained by following the scheme presented 
above. The theory can be seen to follow from imposing the constraints as conditions 
on quantum states, by solving the constraint (101) but keeping the constraint (102). 
Alternatively it can be found by imposing the temporal gauge A = 0, see e.g. [56, p. 
398]. 

In his pilot-wave approach Valentini introduced a field for which the dynamics is 
determined by the guidance equation 

A,-" (no) 



The velocity field is the one that enters in the continuity equation for |\& 



2. 



, d 6 x— \W-=-r =0. (Ill) 
dt J 5Ai V SAiJ y ' 

However, instead of regarding the field A^ as the beable, Valentini suggests to regard 
what he calls a "Faraday geometry" as the actual beable, with a Faraday geometry 
being an element of what he calls "electrodynamic superspace" [23, p. 73]. A "Fara- 
day geometry" is hereby understood as an equivalence class of fields connected by the 
gauge transformations (109) and the "electrodynamic superspace" is the space of such 
equivalence classes. This is analogous to Wheeler's notion of "superspace" in geometro- 
dynamics [85, 86], which is the space of 3-geometries, a 3-geometry being an equivalence 
class of 3-metrics connected by spatial diffeomorphisms. 

Let us consider the ontology in more detail. First note that that the constraint (108) 
is solved by wave functionals that only depend on the transverse part of the vector 
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potential and that such wave functionals obey the functional Schrodinger equation (95). 
As such, Valentini's guidance equation (110) reduces to 




5S 




5S 



= 0. 



(112) 



5M 



Because the wave functional ^(Af, t) satisfies the same functional Schrodinger equation 
as in Bohm's pilot-wave model, the dynamics of the transverse field is the same as in 
Bohm's approach. The longitudinal field is static and has no effect on the dynamics of 
the transverse field. Hence, the longitudinal field appears to be an unphysical degree of 
freedom. Valentini's suggestion to take the equivalence class of fields A; t as the actual 
beable in effect eliminates this unphysical degree of freedom by treating all its possi- 
ble values as physically equivalent. 12 As such, it is clear that, while the ontology in 
Valentini's approach is strictly speaking different from that in Bohm's approach, they 
are empirically equivalent, even in quantum non-equilibrium. 

Of course there exist other ways to eliminate the gauge symmetry (109). For example, 
one could choose a unique representant from each equivalence class, or one could choose 
an invariant function of the field A^ such as the field Bi = €ijkdjA k . While all these 
choices correspond to different ontologies, they would lead to empirically equivalent 
theories, even in quantum non-equilibrium. Such a type of ambiguity does not only 
appear for pilot-wave theories, but is common for gauge theories. For example, in 
classical electromagnetism different possible ontologies are possible that are empirically 
equivalent, such as for example an ontology given by the electric and magnetic field or 
an ontology given by gauge-fixed potential (see [87] for an extensive discussion). 

We finish this section with some notes on the quantum equilibrium distribution in 
Valentini's approach. From (112), it is clear that an equivalence class will evolve to 
an equivalence class. A measure on the set of transverse fields naturally determines 
a measure on the set of equivalence classes. Hence the quantum equilibrium measure 
\ty(Aj,t)\ 2 VA T in Bohm's pilot -wave approach naturally defines a quantum equilib- 
rium measure on the set of equivalence classes. More generally, one could define the 
equilibrium measure using the Faddeev-Popov procedure. Given a suitable gauge fix- 
ing f(Ai) = 0, one has the measure \^(A i ,t)\ 2 5(f(A i ))A F p : f(A, i ) Ylj^Aj on the space 
of fields Ai, which defines a natural measure on the space of equivalence classes. The 
resulting measure is independent of the choice of gauge fixing, so that, in particular, it 
agrees with the one defined in terms of transverse fields. 

In the following sections we will discuss scalar quantum electrodynamics and the 
Abelian Higgs model. We will present a pilot-wave model for these field theories by 
introducing beables only for gauge independent degrees of freedom, thereby avoiding the 
surplus structure in Valentini's approach. For Valentini's approach to scalar quantum 

12 Note that the dynamics of the field A, (including the field Af) is uniquely determined by the 
equations of motion and the initial conditions. Therefore, strictly speaking, this pilot-wave theory does 
not constitute a gauge theory in the sense described in section 2.1.1. Hence, although the transformation 
Ai — > Ai + d{9 is a symmetry of the theory, it is is not a gauge symmetry in this sense. Nevertheless, 
as discussed above Af is clearly an unphysical degree of freedom and can be dismissed on this basis. 
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electrodynamics, see [23]. For his approach to the Abelian Higgs model, as well as 
Yang-Mills theories, see [25]. It is expected that the equivalence between Valentini's 
approach and an approach in terms of gauge independent degrees of freedom, which was 
established here in the case of the free electromagnetic field, can be established much 
more generally. 



5 Scalar quantum electrodynamics 

In this section we present a pilot-wave model for scalar quantum electrodynamics (a 
quantized scalar field interacting with a quantized electromagnetic field). The model 
is obtained by introducing beables for gauge independent degrees of freedom, just as 
in Bohm's pilot-wave model for the free electromagnetic field (for Valentini's pilot-wave 
approach to scalar quantum electrodynamics, see [23]). 

The Lagrangian and the corresponding equations of motion for classical scalar elec- 
trodynamics are given by 

L = Jd 3 x {[D^YD^ - m 2 0> - V(<f>*<f>) - \f^F^ , (113) 

D M Z^ + m 2 + |^ = O, jD ^>*+m 2 0* + ^ = O, (114) 

dpF** = s v = ie {<j>*D v (j> - <pD v *<f) , (115) 

with _D M = dfj, + ieA^ the covariant derivative and s M the charge current. The theory has 
a local U(l) symmetry 

->■ e ief? , 0* -> e~ iee 0* , A" -> N 1 - d^6 , (116) 

which signals the presence of gauge freedom. 13 Hence, just as in the case of the free 
electromagnetic field there will be first class constraints in the Hamiltonian picture. 
Indeed, in the Hamiltonian picture we have [43, pp. 113-127]: 

SL SL 

n^ = - J = D; ( j ) *, n r = -^ = J D o 0, (ii7) 

0(f) 0(fi* 

SL SL 

n Ao = -^ = o, n Ai = - T = (A l + d t A ), (us) 

0A dAi 



13 Unlike the case of free electromagnetism, not all of the transformations (116) form gauge trans- 
formations in the sense defined in section 2.1.1. Writing g = e le0 the local U(l) transformation (116) 
becomes <j> — > g<j>, <j>* g^ 1 ^*, — > + ig^d^g/e. Assuming as boundary conditions that the 
fields A^ vanish sufficiently fast at spatial infinity, the only functions g that preserve these boundary 
conditions are those that have a limit in U(l) at spatial infinity and reach that limit sufficiently fast. 
One can easily show that, of those transformations, only those that go to one at spatial infinity are 
gauge transformations [79]. Modding out over the gauge transformations yields a residual global U(l) 
symmetry. This residual symmetry is then present when the theory is expressed in terms of gauge 
independent variables. 
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h = J d 3 x (jvn^ + (A</>)* A0 + m 2 <f>*(f> + v{4>*4>) 

+ -n Ai n Ai + -F^Fij + a (diU Ai + So ) J , (H9) 
xi = u Ao , X 2 = diU At + s = ftiu. + ie (0*n^* - 0ig , (120) 

where the constraints Xi an d X2 are first class. 

Since U^ is a constraint, the variable A can, just as in the case of the free electro- 
magnetic field, immediately be identified as a gauge variable. In order to identify the 
remaining gauge variables we first perform the following canonical transformation 

= e ie ^0, n^ = e- ie v^%, 

A = A, u At = U At -^s , (121) 

where s"o = ie (0*11^* — ^IT^) is the charge density expressed in terms of the new vari- 
ables. The Hamiltonian for the new variables is given by 14 

H = J d 3 x(u^ + (Dj4>)*Dj4> + m 2 0*0 + V{4>*4>) - \s ^s 

+ ^U Ai U Ai + l -F l3 F l3 + Aoc^n^ , (122) 

where 

Df = $ - ie [b l3 - A 3 =d i - ieAJ (123) 

only depends on the transverse degrees of freedom of A\ and F i3 = diA 3 — d 3 Ai. The 
term containing sq is the Coulomb energy. 

In terms of the new canonical variables the remaining constraint X2 = 9*11^ +sq = 
reads dJl At = 0. So the remaining constraint now has the same form as in the case of 
the free electromagnetic field. By performing a Fourier transformation on A\ and IT^., 
similarly as in the case of the free electromagnetic field, cf. (88), the remaining gauge 
degree of freedom can be separated from the gauge independent degrees of freedom. The 
quantization then proceeds in the same way as in the case of the free electromagnetic 
field. 

The fact that the fields Aj, 0, 0* and the corresponding conjugate fields are gauge 
independent variables is well known, see e.g. [81, pp. 302-306], and motivated the canon- 
ical transformation (121). The resulting theory can also be obtained by imposing the 
Coulomb gauge. 



14 Note that, as alluded before in footnote 13, there is a residual global 17(1) symmetry, namely 
—> g4>, <f>* — > <? _1 </>*, — > — > ^n^, , with g a constant phase factor. 
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After quantization and by using the functional Schrodinger picture, we obtain the 
functional Schrodinger equation 

= / Sx ( - w» + ( w Dji + mm + vm 



2 H* T 6(f)J V 2 
1 5 2 1 



26AJAT 2 A'VM;j*, (124) 

for the wave functional \l/(0, 0*, , t). Beables 0, 0*, can be introduced which evolve 
according to the guidance equations 
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= nr (125) 

and for which the quantum equilibrium density is given by |\I/(0, 0*, AJ , t)\ 2 . 



6 Massive spin-1 field 

In this section we present a pilot-wave model for a massive spin-1 field. The reason to do 
this, is that we will encounter a massive spin-1 field when we discuss the Abelian Higgs 
model in the next section. The pilot-wave model was presented before by Takabayasi 
[21]. 

The Lagrangian for a classical spin-1 field with mass m is obtained by adding a mass 
term to the Lagrangian for a classical massless spin-1 field (which is just the Lagrangian 
for the electromagnetic field) 

L — Jd 3 x (--^F^ + . (126) 

The corresponding equation of motion is 

d^ v = -m 2 A v . (127) 
The Hamiltonian formulation is determined by [43, p. 35] 

SL SL 

nA = -^ = o, n Ai = 71 - = (A + 9 i A ) ), (128) 

oA Q dAi 

/( 1 I 777^ 777^ \ 

d 3 x I -U Ai U Ai + -F l3 F tJ + —AiAi - —A\ + A diU Ai \ , (129) 
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Xi = n A , X2 = diU Ai - m 2 A , (130) 

where Xi and X2 are second class constraints. By performing the canonical transforma- 
tion 

A = A + -^n^ , n Ai = u Ai , (131) 

we find that the constraints reduce to A = U Aq = 0. The unconstrained variables are 
A and n^.. The Hamiltonian for these variables reads 

H = jd*x + -L (d^) 2 + iFyFy + ^AA^j . (132) 

After quantization and by using the functional Schrodinger picture, we obtain the 
functional Schrodinger equation 

i? = / <? X {- l -J^ - J_ ( a i *) 2 + ^F l3 F l3 + —AiAi I * . (133) 
dt J y 2 5AiSAi 2m 2 \ l 8Aj 4 13 13 2 1 ! 

The guidance equation for the field beable A reads 

7 i5S 1 „„ 

5ij m J SAj 

7 Abelian Higgs model 

The electromagnetic interaction is unified with the weak interaction in the electroweak 
theory of Glashow, Salam and Weinberg. Instead of introducing beables for the elec- 
tromagnetic field, it might therefore be more fundamental to introduce beables for the 
electroweak field. However, the electroweak theory is a non- Abelian gauge theory, based 
on the symmetry group 577(2) x C/(l)y. Due to the technical difficulties that come with 
the construction of a pilot-wave model for non- Abelian gauge theories (cf. section 8), we 
do not make an attempt to construct a pilot-wave model for the electroweak theory. 

However, there is one important aspect of the electroweak theory that we want to 
discuss in the context of pilot-wave theory, namely spontaneous symmetry breaking 
and the associated Higgs mechanism. In the electroweak theory, the SU{2) x U(l)y 
symmetry gets broken to the U(l) em symmetry of electromagnetism. In this process the 
weak interaction bosons W + , W~ ', Z° acquire mass through the Higgs mechanism. That 
the pilot-wave theory has no problem in dealing with spontaneous symmetry breaking 
and the Higgs mechanism can be illustrated in the simple case of the Higgs model (which 
has a U(l) symmetry) (for Valentini's pilot-wave approach see [25]). 

Before turning to quantum field theory, let us first consider the Abelian Higgs model 
on the level of classical field theory. The Abelian Higgs model (see e.g. [88]) describes 



42 



a scalar field coupled to an electromagnetic field, for which the classical equations of 
motion can be derived from the Lagrangian 

L — Jd 3 x {{D^YDy + ^ 2 0*0 - A(0*0) 2 - \f^F^ . (135) 

The constant /j, is real, so that the mass of the scalar field is imaginary, making the field 
tachyonic. The constant A is real and positive. This Lagrangian is the one of scalar 
electrodynamics, cf. section 5, with potential V{i^)*<p) = A(0*0) 2 . The theory exhibits a 
local U(l) symmetry, given by the transformations in (116). 

Broadly speaking, one has spontaneous symmetry breaking if the ground state (the 
state of lowest energy) does not share the symmetry of theory, but is degenerate in such 
a way that the different ground states are mapped to each other by the symmetry. To 
each ground state there corresponds a set of solutions to the equations of motion that 
arise as perturbations around that ground state. The effective equations of motion for 
each of those sets generically do not have the symmetry of the original equations of 
motion. 

With the spontaneous symmetry breaking of a gauge symmetry new features arise, 
such as the generation of mass through the Higgs mechanism. However, since gauge 
symmetries are unphysical, some care is required with the meaning of symmetry breaking 
in this case. This will be considered in detail in [79]. Strictly speaking, the symmetry 
that is broken in the case of the Abelian Higgs model is not the gauge symmetry but the 
residual global U(l) symmetry that arises after modding out the local U(l) symmetries 
by the actual gauge symmetries (see footnote 13). 

One way to see this is by considering the formulation of the theory in terms of gauge 
invariant variables for which the Hamiltonian is given by 15 

H = J d 3 x(u r U^ + (Dj4>)* Dj4> - /i 2 0*0 + A(0*0) 2 - \s Q ^s Q 

+ \u a t U Af - \AjV 2 A^j , (136) 

where we dropped the bars on top of the fields for notational simplicity. 

In this theory the gauge symmetry has been eliminated, but, as noted before, it 
possesses a residual global U{1) symmetry 

0->e ief V, 0*->e" ie V, U^e' ie9 U^, U r -> e ied U r . (137) 

The minima of the energy, i.e. the Hamiltonian, are given by = a/ /i 2 /2Ae 1Q , 0* = 
a/ /i 2 /2Ae~ ia (where a is constant) and = n^. = Aj = U a t = 0. Let us now choose 
one of those ground states, say = a/ /i 2 /2A = v / \/2, and make the following change of 

15 Usually the analysis of symmetry breaking in the Abelian Higgs mechanism is carried out by 
employing the unitary gauge. 
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variables (which forms a canonical transformation) 



<t>=fyv + ri + %), = -j={n v - in 5 ), 

0* = -L( u + ?7 -^), n r = -L(n, + in f ), (138) 

where 77, 11^, £, 11^ are real variables. Let us also introduce 

Af = —d£, II a l = ev^U^. (139) 



ev 



This transformation is not canonical; the Poisson bracket of A\ and U a l is given by the 
longitudinal delta function. The latter variables allow us to introduce the variables 

A t = Aj + Af, Tl At = Tl A T + Tl A L. (140) 

The variables A^ \i A . form canonical pairs again. 

Expressing the Hamiltonian in terms of the variables 77, U v , £, LT^, A iy H Ai , only up to 
quadratic terms and up to a (infinite) constant, we obtain 

+ \^ Ai + ^ + —^-AiAi + -^F^ . (141) 

The higher order terms contain the coupling between the fields. The first line is recog- 
nized as the Hamiltonian of a scalar field with mass a/2/j 2 and the second line as that 
of a spin-1 field with mass Ve 2 v 2 . What has happened is that the Goldstone boson 
£ and the transverse field Af together form a massive spin-1 field. This is the Higgs 
mechanism. 

So far the discussion of the symmetry breaking and the Higgs mechanism was done 
on the classical level. A similar picture arises in the pilot-wave approach to the cor- 
responding quantum field theory. Since the Lagrangian of the Abelian Higgs model 
corresponds to the one of scalar quantum electrodynamics, the quantized theory and 
the corresponding pilot- wave model are those presented in section 5. The functional 
Schrodinger equation reads 



I d3x {'W^ + ( W ^ - + A (<^) 2 

^-4)^(^-4)-^-^)^ (l42) 
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and the guidance equations are 



• 5S 2± 1 fJS 5S 
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(143) 



Let us now use the expansion = (v + r\ + i£)/y/2 (just as in (138)) in the functional 
Schrodinger equation and in the guidance equations. We will write the functional Schro- 
dinger equation only up to quadratic terms in the fields r], £ and Af, and up to a (infinite) 
constant. The guidance equations will only be written up to linear terms in the fields. 
The higher order terms represent the interactions between the fields. By introducing the 
field Ai = Aj + Af, where we define Af = —d^/ev, and by making use of the relations 

5 1/5 5\ 5 5 5 5 d t 5 , AA , 

+ T77:> T7T = ^-4-, (144) 



50 >/2\^ ' 5A i 5AJ 5Af SAf V 2 5£' 



we find 



M f j3 / 1 5 2 1 . . 2 2 1 

i^TT = / « - -— + -Oi-qdiT] + (J, T) 



dt J V 2 5r/ 2 2 " * ' ' ^ ' 2 5A i 5A i 



2eV 

Hence we recognize a scalar field 77 with mass a/2/x 2 and a spin one field with mass 
V e 2 v 2 . The guidance equations reduce to 

• SS ■ 5S 1 00 55 

" = ■*;• A '-iA,-^ did 'JA- (146) 

which are the familiar guidance equations for a real scalar field and a massive spin-1 
field. 



8 A note on non-Abelian gauge theories 

For the development of a pilot-wave model for a gauge theory, our approach is to intro- 
duce only beables for gauge independent degrees of freedom. As explained before, the 
way to do this is by starting from the formulation of the quantum field theory which 
is obtained by quantizing only gauge independent degrees of freedom. For the gauge 
theories we have considered so far, this approach worked very well. One reason for this 
is that it was easy to identify and isolate the gauge independent degrees of freedom as 
canonical pairs. For non-Abelian Yang-Mills theories this is much harder. The reason is 
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that the constraints are non-linear in the fields. For example for the free Yang-Mills field 
the constraints read n^g = and D^Ha? = c^H^a — f abc A^IlA<r = 0. These constraints 
look similar to the constraints for the free electromagnetic field, with the difference that 
the second constraint contains a non-linear term. Some progress has been made on the 
identification on gauge independent degrees of freedom for non-Abelian gauge fields, see 
e.g. [89] and references therein. This work could potentially be used for the development 
of a pilot-wave model. 

In Valentini's approach [25] , where the actual beable is an equivalence class of gauge 
equivalent fields, no difficulties appear on the level of the dynamics, but instead on the 
level of specifying a quantum equilibrium distribution. In specifying the equilibrium 
distribution distribution, one needs to avoid the difficulties associated with the infinite 
gauge volume. As in the case of the electromagnetic field (cf. section 4.3), this could 
be done by applying the Faddeev-Popov procedure. However, it was shown that a large 
class of gauge fixings are ruled out [83, 84]). 

9 Fermionic fields 

In this section we consider two pilot-wave approaches for fermionic fields that introduce 
fields as beables. One approach is by Holland [17, pp. 449-457] and [22], the other by 
Valentini [23,24]. We point out some problems with both these approaches. 

9.1 Holland's approach: Euler angles as beables 

Before turning to Holland's pilot-wave model for fermionic fields, we first have a look 
at a particular pilot-wave model for non-relativistic spin-1/2 particles, which was also 
proposed by Holland, cf. [17, pp. 424-448] and [22,90]. It is instructive to first con- 
sider this model because it bears some resemblance to his model for fermionic fields. 
The formulation of both models involves angular variables, but these variables play a 
significantly different role in each of the models. 

9.1.1 Non-relativistic spin-1/2 particles 

The standard wave equation for non-relativistic spin-1/2 particles is the Pauli equation. 
Holland's starting point is a different wave equation which would underly the Pauli 
equation (although there is no evidence for such an underlying theory for elementary 
particles). 

For simplicity, we restrict our attention to a single spin-1/2 particle. In this case, the 
quantum state is given by a wave function ck, t), defined on the configuration space 
R 3 x 577(2). As before, the coordinates x represent the position degree of freedom. The 
coordinates a = (a,/?, 7), with 0^a^7r, 0^/3< 2n, ^ 7 < 47r, are Euler angles 
that parametrize SU (2) and represent the spin degree of freedom. 

The group SU (2) is introduced as the double covering group of the Euclidean rota- 
tion group SO (3). The angles ex. stand in two-to-one relation to a set of Euler angles 
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parametrizing 5*0(3), denoted by a' = (a', f3', 7'), with ^ 0/ ^ 71, ^ j3' < 2n, 
^ 7' < 2n, through the relations a' = a, (3' = (3, 7' = 7 mod 2n. 

The Schrodinger equation for this wave function is obtained by quantizing the classi- 
cal rigid rotator (for which the configuration space is 1R 3 x SO (3)). The Pauli equation, 
which is the standard equation describing non-relativistic spin-1/2 particles, can be ob- 
tained from it. Actually, because of the use of the double covering group, this theory 
contains the Pauli equation twice [17, pp. 439-441]. 

The beables in Holland's model are a position x and an orientation a' e SO (3) 
in 3-space. The orientation a' is obtained from a beable configuration a e SU(2). 
The position x represents the center of mass of a rigid homogeneous sphere and a' 
represents the orientation of a fixed axis through the center of this sphere. The exact 
form of the Schrodinger equation and the guidance equations are not really important 
here. The equilibrium distribution for the beables (x, a.) is given by |^(x, et,t)\ 2 d 3 xdQ, 
where dVl = sin adadfid^ is the measure on SU(2). 

This pilot-wave model has the property that any two wave functions are always very 
much overlapping in a-space. This is for example illustrated by considering states of 
the form: 

■0(x, ot,t) = i/j + (-x.,t)u + (ac) + ^-(x, t)u-(ot) , (147) 

where 

u + = — ^- cos(a/2)e- i(/3+7)/2 , m_ = -^L- sin(a/2)e i(/3 - 7)/2 . (148) 
2V2n 2y/2n 

The states u + and w_ respectively correspond to spin up and spin down along the x%- 
direction and forms the Pauli spinor. States of this form already contain the 

full Pauli theory. The fact that different wave functions always have much overlap in 
cn-space implies that the beable (X, or a', never records any outcome of a measurement, 
in particular it does not record the spin of the system. Of course we are considering only 
a single particle here and there is no need for measurement results to be recorded in the 
beables of one particle. On the other hand, if we perform a spin measurement with a 
Stern-Gerlach device, then the spin will get recorded in the position beable x, since the 
wave function will develop non-overlapping terms in position space. For example, if we 
perform a spin measurement in the ^-direction, then the wave functions ip+ix^t) and 
ijj- (x, t) will separate in the a^-direction. As such, the position beable x will either move 
up or down in the a^-direction, thereby displaying the outcome of the spin measurement. 

Considering the many particle case, it is unclear to what extent measurement results 
get recorded in the configuration of angular beables. However, this is not really impor- 
tant if the goal is merely to show that the pilot-wave model reproduces the standard 
quantum predictions, since measurement results will at least be recorded in the positions 
of the rigid spheres. 

The insight that measurement results generally get recorded in positions of things, 
like the positions of measurement needles, led Bell [91, 92] to formulate a more natu- 
ral pilot- wave model for non-relativistic spin-1/2 particles where the only beables are 
particle positions (starting from the Pauli equation). 
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9.1.2 Fermionic Schrodinger field 

Let us now turn to Holland's pilot-wave model for fermionic fields. Holland illustrated 
the construction of his model for the case of the fermionicalfy quantized Schrodinger 
field, but his method applies equally well to other fermionic quantum fields, like for 
example the quantized Dirac field. 

Functional Schrodinger picture The fermionically quantized Schrodinger field is 
obtained by imposing the anti-commutation relations 

Wx),f(y)} = *(x-y), {#(x),#(y)} = W>(x),V>(y)} = 0. (149) 

The corresponding Hamiltonian operator has the same form as the bosonic Hamiltonian, 
which was given in equation (54) (the term E was not present in Holland's presentation). 

Holland chooses a representation for these field operators in terms of angular vari- 
ables («k,/3k,7k) in momentum space. 16 For each momentum k, the angular variables, 
collectively denoted as a k = (a k , /?k, 7k), with range ^ a k ^ n, ^ /3k < 2n, 
^ 7k < 47r, parametrize the group SU(2). Holland further assumes a finite spatial 
volume so that the set of momenta is discrete. 

In this representation the field operators act on wave functionals ty(at,t). The wave 
functionals can be written as superpositions of states of the form Yl k Wkn k (c*k), where 
each Mkn k (cKk), with = ±, is given by (148) for each momentum k. For example, the 
vacuum state is given by 

«k-(«k) ■ (150) 

k 

States describing n particles with definite momenta k 1; . . . , k n are given by 

*k!...k„(Q:) = s(ki,...,k n )u kl+ (a kl ) . . . ii kn+ (a k J JJ « k -(a k ), (151) 

k^ki,...,k„ 

where s(k 1; . . . , k n ) is the sign of the permutation that maps (k l5 . . . , k n ) to (k' 1; . . . , kJJ, 
where {k 1; . . . , k n } = {k' 1; . . . , k' n } and k' : < • • • < k' n (see footnote 16 for more on 
the order relation on the set of momenta). (More precisely we have \I/ki...k n (aO = 
(o:|a kl • • .a^l^o), where a k are the fermionic creation operators of particles with mo- 
mentum k. Using [17, pp. 449-457] and [22], it can easily be verified that this expression 
reduces to the one in (151).) 

16 We do not reproduce the equations of Holland's model here. There is however one small correction 
we want to suggest to Holland's presentation. In order to express the fermionic creation and annihilation 
operators a k and ak in terms of angular momentum operators, a total order relation ^ is needed on 
the set of possible momenta [93]. This order relation can be chosen completely arbitrary. However, the 
order relation Holland specifies is only a partial one and not a total one, i.e. not every two momenta 
can be compared according to Holland's order relation. Holland's order relation is defined by: k' < k 
if k\ < ki for one i — 1,2,3, and k'- < kj for j ^ i, which does not yield a total order relation. However 
a total order relation can easily be found. One could for example consider the one defined by: k' < k 
if k[ < ki for one i — 1,2,3, and kj ^ kj for j > i. However, this correction does not imply any 
formal change to Holland's presentation. It should just be kept in mind that some particular total 
order relation should be chosen. 
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Ontology Holland starts with introducing beables for the operators ct^. The 
quantum equilibrium distribution for the beables is given by [^(a;)) 2 Y[ k dVtk- However, 
these beables form a field in momentum space and unfortunately Holland does not 
explicitly specify how these beables relate to some ontology in space. 

When discussing the classical limit of his model [17, pp. 455-457], Holland defines 
a complex-valued field ijj(x.,t) in terms of the beables and the wave functional, for 
the special case of a vanishing quantum potential (in [17] Holland adopts, just as Bohm 
for the electromagnetic field, a second order view in terms of the quantum potential). 
Possibly a similar definition can be given in general (i.e. also when the quantum potential 
does not vanish). One way to do this could be by using Holland's local expectation value 
for the operator ^(x). 17 

Another possibility would be to introduce a beable for the particle number density 
operator ■0^(x)'0(x), again by means of Holland's local expectation value. 

Reproducing the quantum predictions Whether or not Holland's model is capa- 
ble of reproducing the standard quantum predictions depends of course on the particular 
choice of ontology in physical space. Macroscopically different states should correspond 
to different beable configurations. Therefore it is necessary (at least for the type of be- 
able suggested above) that the wave functionals are non-overlapping in the configuration 
space of angular beables. This guarantees that macroscopically different states corre- 
spond to different beable configurations a. It is then further necessary that the different 
beable configurations a. correspond to different beables in physical space. Since our ex- 
amples of possibly ontologies in physical space are rather complicated constructions 
involving the wave functional ^(ck), it is not immediately clear that this is guaranteed. 

It is easy to see that states containing only a small number of particles have much 
overlap. However some states have little overlap. Consider for example the state 
^ki...^^)- This state will have little overlap with the vacuum state \E'o( Q: ) i n the 
case the number of particles n is high enough. This follows from the fact that the 
overlap of the functions Wki+^ki) • • • Wk n +(o: kn ) and Wki-(£*ki) • • • Wk n -(o:k n ) decreases 
when n increases (this can be seen by considering the law of large numbers 18 ). Similarly, 
two states ^^...^(ot) and ^^...k' , (<*) will have negligible overlap when the number of 
momenta in the symmetric difference of the sets {ki, . . . , k n } and {ki, . . . , k n } is suf- 
ficiently large. It is an open question whether states describing localized macroscopic 
systems have negligible overlap. 

So it is as yet unclear whether Holland's model is capable of reproducing the pre- 
dictions of standard quantum theory. In the following section an alternative but similar 
model is considered for which this question is much easier to analyse. However, the 

17 Holland only introduced local expectation values for Hcrmitian operators [17]. One could define 
the local expectation value for ip(x) as the linear combination of the local expectation values for the 
Hermitian operators Vv(x) = (V>(x) + ^>t(x))/2 and ^»(x) = — i(^(x) — ^>t(x))/2. Note that the local 
expectation value for V>(x) would be a complex- valued field and hence would not satisfy the anti- 
commutation relations (149) of the field operator V>(x). 

18 I am grateful to Owen Maroney for this observation. 
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upshot of that analysis is that the model seems inadequate. 

9.1.3 Alternative model for the fermionic Schrodinger field 

Since it is rather hard to analyse the viability of Holland's model for the fermionic 
field due to its complicated ontology, it is worthwhile to consider the following simpler 
alternative. In developing his pilot-wave model, Holland started with the field operators 
ak and a k in momentum space, which satisfy the anti-commutation relations {a^, a k ,} = 
<5kk'- By using a representation in terms of angular momentum operators in momentum 
space, he was led to introduce the beables a?k- Instead of starting with the field operators 
in momentum space, one could also start with considering operators in physical 3-space. 
Assuming a spatial lattice, instead of a continuum space, the anti-commutation relations 
for the field operators ^ x and ^* read {-?/> x , ■?/>*,} = 5 XX /, where x and x' are lattice 
sites. These anti-commutation relations are formally the same as the anti-commutation 
relations for the field operators in momentum space. 19 One could therefore adopt the 
techniques introduced by Holland and use a representation of these ant i- commutation 
relations in terms of angular momentum operators. This time the angular momentum 
operators would be operators in physical space. This would result in beables ai x = 
(a x ,/3 x ,7 x ) in physical space. While one could further assume a rigid body at each 
lattice point x, whose rotational motion is described by the angles o: x (similarly as in 
Holland's model for non-relativistic spin-1/2 particles), there is no need for that. 

We do not present the dynamics of the beables here. It could easily be found by using 
the results in [47] . Even without knowing the dynamics explicitly we can say something 
about the viability of the model. This could be done by comparing the typical beable 
configurations corresponding to a state \& m which represents matter in a certain region 
V in physical space and a state \l/ e which represents an empty region V. We have, in 
analogy with (150) and (151), that 

* e (a) = n «x_(a x ) (152) 

x in V 

and 

*jif(a) = s(xi, . . . ,x n )u xi+ (a xi ) . . .u Xn+ (a x J J | w x _(a x ), (153) 

x in V 
x^xi ,...,x n 

where x 1; . . . , x n are the positions of the particles in V. For our purposes the environ- 
ment can be ignored in the states. 

Note that matter will never be distinguished from empty space from the beables f3 x 
or 7 X , because the states states ^ m and \l/ e yield the same distribution for them. They 
could only be distinguished from the beables a x . This could be done by considering 
the quantity Ay = ^ x in v a x / n. If this quantity is typically different for the states 

19 Holland carried out his analysis for a discrete set of momenta. However, probably a similar (for- 
mal) analysis can be carried out in the case of a continuous set of momenta. In that case we could, 
equivalently, assume a physical space which is continuous instead of discrete. 
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m and \I/ e then one could read off from the beables whether the region V is filled 
with matter or not. For a state \& the expected value for the quantity Ay is given by 
(A v )q, = J A V \^J\ 2 rj x in v <ifi x . Matter can be distinguished from empty space from the 
beables if 

\{A v )v u - {A v )v e \ > max(A* M A y , A* e A v ) , (154) 

where Aq,A v = a/ (A v )q, — (A v )% is the standard deviation of Ay for the state For 
a cubic region V with edges of length L, matter with particle number density p and a 
lattice spacing a, we find in Appendix B that (154) reduces to 

Lap 2/3 > 1 . (155) 

Assume one particle per atomic distance cubed, so that p = 10 30 /m 3 . Then (155) 
becomes La ^> 10~ 20 m 2 . If we further assume a lattice spacing of the order of the 
Planck length, i.e. a = lCT 35 m, then L 3> 10 15 m. This means that in this case we can 
only distinguish matter from empty space by considering length scales that are much 
bigger then 10 15 m. This is clearly unreasonable. The lattice spacing could of course be 
taken much larger than the Planck length in a non-relativistic theory. But presumably 
it should not be taken much higher then 10~ 15 m which is only three orders lower than 
the Compton wavelength of the electron. For the latter value of the lattice spacing, 
L must be much bigger then 10 _5 m. This is still a rather large distance, as we can 
easily distinguish such length scales with the naked eye (the diameter of a human hair 
typically ranges from 10~ 5 m to 10~ 4 m). This implies that it seems very unlikely that 
the beables in this model yield an image of general macroscopic matter distributions. 
As such the model seems to fail to reproduce the standard quantum predictions. 

There also seems to be little hope that an application of this model to relativistic 
quantum field theory will yield a better result. In a relativistic model the beables would 
consist of four set of angles at each point in space (one corresponding to each component 
of the Dirac spinor) while the lattice spacing should be taken much smaller. 

9.2 Valentini's approach: Grassmann fields as beables 

In [23, 24] Valentini suggested a pilot- wave model for fermionic fields in which the beables 
are anti-commuting fields, also called Grassmann fields. In this section we point out 
a problem with this approach. We will see that although it is possible to introduce 
dynamics for anti-commuting fields, it is unclear what the actual guidance equations 
and corresponding equilibrium distribution should be, in order to have a pilot-wave 
model that reproduces the standard quantum predictions. 

The starting point of Valentini's approach is to write the fermionic field theory in a 
functional Schrodinger picture in which the quantum states are functionals defined on 
a space of anti-commuting fields. Although this functional Schrodinger picture, which 
was initiated by Floreanini and Jackiw [94], bears some resemblance to the functional 
Schrodinger picture for bosonic field theories, there are some important differences. 
One is that the fermionic wave functionals are not complex-valued, but take values 
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in a Grassmann algebra. A consequence is that the wave functional does not yield a 
transition amplitude for finding the state in a certain field configuration. 

In order to point out clearly the problem with Valentini's model, we present the 
functional Schrodinger picture for fermionic fields in some detail in section 9.2.1. The 
reader who is familiar with this picture can immediately move to section 9.2.2. 

As mentioned in section 3.3, Valentini presented his pilot-wave model for the quan- 
tized Van der Waerden field. Valentini used the Van der Waerden theory, instead of the 
equivalent Dirac theory, in order to circumvent the problem of the constraints in the 
Dirac theory. For simplicity we do not discuss Valentini's pilot-wave model for the quan- 
tized Van der Waerden field, but for the fermionically quantized Schrodinger field (which 
was already introduced in section 9.1.2). The formulation of the functional Schrodinger 
picture for the quantized Schrodinger field is completely similar to that for the quantized 
Van der Waerden field or the quantized Dirac field. The functional Schrodinger picture 
for the Dirac field can for example be found in [95] . 



9.2.1 Functional Schrodinger picture in terms of Grassmann fields 

In order to set up the functional Schrodinger picture for the fermionically quantized 
Schrodinger field, we first introduce some basic elements of Grassmann algebras and 
analysis over these algebras. Details can be found in [96-98]. 
Consider two fields u>(x) and iu^(x) which satisfy 

{w(x),w(y)} = Mx)V(y)} = Mx),u%)} = 0. (156) 

These fields are Grassmann fields (i.e. for each x the values tu(x) and uA(x) are Grass- 
mann numbers). The algebra generated by these fields forms an infinite dimensional 
Grassmann algebra, which we denote by A. Its elements are functionals G(w, w^) of the 
form 



G = G 00 + J rf 3 xiGio(xi)w(xi) + J d 3 x 1 G i{xi)w t (x 1 ) + ... 

+ J d 3 x 1 . . . d 3 x n+m G nm (x 1 , x„ +m )w(x!) . . . w(x n )w t (x„ +1 ) . . . w f (x n+m ) 
+ . . . , (157) 

where Goo is a complex number and where G nm (xi, . . . ,x n+m ), with n + to > 0, are 
completely anti-symmetric complex-valued distributions. There is an involutive map 
G i— ?■ G^ on A, called complex conjugation, which is linear under addition, anti-linear 
under multiplication by complex numbers, satisfies (GiG 2 Y = G\G\ and (u>(x))* = 
u>t(x), and which acts as the ordinary complex conjugation on complex numbers. An 
element G of A is called odd if G nm = for all (to, n) for which n + to is even. All odd 
elements anti-commute among themselves. The subspace of odd elements is denoted by 
A„. 
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We can now introduce the space J 7 of A- valued functionals F(u, w), which are defined 
on the space of A„- valued fields -u(x) and ■u'l'(x), 20 and which are of the form 21 

F = F 00 + J rf 3 a; 1 F 10 (x 1 )M(x 1 ) + J d 3 x 1 F 01 (x 1 )^(x 1 ) + . . . 

+ J d 3 x 1 . . . d 3 x n+m F nm (x 1 , . . . , x n+m )w(xi) . . . m(x„)m 1 '(x„ + i) . . 

+ ..., (158) 

where F 00 is a complex number and where F nm (x 1; . . . , x n+m ), with n + m > 0, are 
completely anti-symmetric complex-valued distributions. 

For functionals in J 7 the functional derivatives 5F(u, u^)/5u(x.) and 5F(u,u^)/5u(x) 
exist for all A G -valued fields it(x) and ■u'(x). The functional derivative 8F(u,w)/5u(yi) 
is hereby defined as the coefficient that arises in the variation SF(u, w) under an in- 
finitesimal variation 5u(x) of -u(x) in the space of A„- valued fields [98], i.e. 

6F(u,v!) = [ d 3 x5u(x) SF } U : V ; ) . (159) 
J du(x.) 

The functional derivative SF(u, u')/(w(x) is defined in a similar way. These definitions 
imply that the functional derivatives anti-commute among each other and that they 
satisfy the properties 

«(x), = l«t( x ), -4--} = 6(x - y) , (160) 



' My) J I ' $ u Ky) 

" (x) -*4)}^W*4)}^°' (161) 

For functionals in J 7 one can also define functional integrals. First of all, for a fixed 
point x, the integral over w(x) is defined by the relations J du(x) = and J du(x.)u(x) = 
1. The integration is hence defined formally, i.e. measure-theoretical concepts play 
no role. The functional integral J VuF(u, v)) is then defined as the (formal) product 
integral J (Y[ x du(x))F(u, v)). The functional integral J Vu^F(u, v)) is defined similarly. 

In the functional Schrodinger picture for the quantized Schrodinger field, the field 
operators ^(x) and V>*(x) are realized by the representation 22 

ftx) -> -L (»(x) + J^) , ?(x) -> ± («t(x) + ^) , (162) 



20 As is standard practice, we treat u(x) and «t(x) as independent fields. Alternatively one could 
consider two real fields Ui(x) and «2(x) (the fact that they are real means that Uj(x) = w|(x), « = 1,2), 
and define u(x) = («i(x) + iM 2 (x))/v / 2 and u^(x) = (ui(x) — iu2(x))/v / 2- 

21 One can also define A Q -valued fields that are not of the form (158) by considering the F nm to be 
A- valued instead of complex- valued. However such functionals are of no interest here. 

22 Just as in bosonic field theory, there exist various other possible representations of the field opera- 
tors. Another possible representation is V>(x) — > u(x), ip*(x) — > S/Su^(x), see e.g. [57, pp. 217-223] or 
[99]. 
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which act on quantum states which are elements of J- ' . The associated inner product is 
defined as 

(^|$ 2 ) = JWDm^ 2 , (163) 
where ty* is the dual of \1/ given by 

^*(u,u j ) = J Vu j Vuexp(uu^ + u j u)¥(u,rf) , (164) 

where we used the notation uirf = J d 3 xu(x)u^(x) and similarly for v)u (the integration 
variables -u(x) and u\x) are A -valued fields in space). The functional Schrodinger 
equation is then obtained by using this representation for the Hamiltonian operator 
(which has the same form as the one for the bosonically quantized Schrodinger field, cf. 
(54)). For the functional Schrodinger representation for the quantized Dirac field one 
just needs to add a spinor index a = 1, . . . , 4 to the Grassmann fields. 



9.2.2 Pilot-wave model 

Let us now turn to the pilot-wave model of Valentini. Valentini postulated anti-com- 
muting, i.e. A G -valued, field beables w(x) and i^(x), which are guided by the wave 
functional ^(w, t) through guidance equations which are of the form u = v^{u,v),t) 
and v) = {y u\ t)Y . The exact form of the guidance equations is not really important 
here. But we do want to note that there are actually some technical problems with those 
guidance equations. They for example contain inverses with respect to the multiplication 
of elements of the Grassmann algebra, but these are not always defined. However, these 
problems, which were discussed in detail in [100], are of minor importance, since, in 
principle, there is no problem for having equations of motion for anti-commuting fields 
(see for example also the generalization of classical mechanics in which some of the 
dynamical variables have anti-commuting values [101-104]). 

The more important problem is that there seems to be no real basis for adopting 
Valentini's guidance equations. Valentini regarded the quantity ^(u, v) , t)^(u, u\ t) as 
the density of fields and found the guidance equations by considering the asociated conti- 
nuity equation. 23 However, this interpretation of seems to be unfounded. The quan- 
tity ^(u, u\ t)^/(u, u\ t) is, just as the wave functional ^(w, u\ t), A-valued. 24 More- 
over, it is unclear what measure should be considered on the space of anti-commuting 
fields. To our knowledge suitable measure-theoretical notions have not yet been intro- 
duced for such a space. Remember also that the integration theory was defined formally 
and was not based on measure-theoretical notions. Therefore seems to have no 
meaning as a density of fields. To be sure, this problem is not a consequence of dealing 

23 Note that Valentini regarded as the density of fields and not ffr*^. Though, regarding 
as the density would be as problematic. 

24 The wave functional ^(u, u\t) should not be confused with a transition amplitude (u, M t |\l/(t)): 
^(u,u\t) is an element of the Grassmann algebra A and (u, u^\^(t)) should be a complex number by 
the definition of the inner product. 
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with a field theory; it would still be present in a regularized theory that describes only 
a finite number of degrees of freedom. So, in summary, it is unclear what the dynamics 
and corresponding equilibrium distribution should be so that the standard quantum 
predictions are reproduced. 

10 A minimalist pilot-wave model 

Motivated by the difficulties of introducing field beables for fermionic field theories, 
Struyve and Westman [26, 27] proposed a model in which there are only beables for the 
bosonic degrees of freedom and no beables for the fermionic degrees of freedom. This 
approach was illustrated for quantum electrodynamics. In that case the beable is a 
transverse field Aj, representing the transverse part of the vector potential, just as in 
Bohm's model for the free electromagnetic field. The dynamics is very similar to that 
of Bohm's model. The difference is that in this case the wave functional also contains 
fermionic degrees of freedom, which have an effect on the dynamics of the beable Aj. 
This makes the beable Aj evolve in such a way as to give the appearance of matter. 

It was argued in detail that this model is actually capable of reproducing the pre- 
dictions of standard quantum field theory, as one can find an image of the visible world, 
and in particular of measurement outcomes, in the field beable Aj. However, while this 
model can arguably be maintained as a logical possibility, it is rather far removed from 
our everyday experience of the world. 

After presenting this radically minimalist model, Struyve and Westman also illus- 
trated how one could introduce additional beables for the fermionic degrees of freedom 
[27]. By means of Holland's local expectation value, beables could be introduced for 
operators such as the mass density operator or charge density operator. In this way 
the beables for the fermionic degrees of freedom are completely determined by the field 
beable Aj and the wave functional (while they do not influence the field beable Aj or 
the wave functional). 

The beables corresponding to these density operators would behave in a way similar 
to that of the mass density in collapse theories [105]. For example, in a double-slit 
experiment with an electron, the mass or charge density of the electron would evolve 
through both slits and would display some interference pattern just before arriving to 
the detection screen, only to "collapse" to an approximately localized density when an 
actual position measurement is performed. This collapse happens vecause during this 
measurement process the wave function of the electron gets entangled with that of a 
macroscopic pointer and with that of the radiation being emitted and scattering off 
the macroscopic pointer. The entanglement with the state of radiation is ultimately 
responsible for the collapse of the mass or charge density of the electron. 
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11 Conclusions and outlook 



In this paper we have studied the construction of pilot-wave models for quantum field 
theory for which the beables are fields. We started with presenting a recipe according 
to which such pilot-wave models could be developed, but only for bosonic fields. For 
gauge theories the recipe implies that beables are only introduced for gauge independent 
degrees of freedom. This recipe can easily be applied in the case of simple quantum 
field theories, as was illustrated for example for the Schrodinger field and for scalar 
quantum electrodynamics. For more involved theories, like Yang-Mills theories, this is 
more complicated and has not yet been achieved. 

In the paper, we did not consider alternative approaches to gauge theories, such as 
for example loop-based approaches. In a loop-based approach the basic variables are not 
the gauge potentials, but loops in space or space-time [106-108]. For example, the free 
electromagnetic field is easily reformulated in terms of such variables [106] . Considering 
that such an approach can also be applied to quantum gravity, namely in the loop 
quantum gravity approach, it might be worthwhile to try to develop a pilot-wave theory 
in terms of such variables. 

For fermionic quantum field theories it is much harder to construct a pilot-wave 
model with fields as beables. We considered a model by Holland and a model by Valen- 
tini. In Holland's model the beables form a field of angular variables. Further study of 
this model is required in order to see whether it reproduces the quantum predictions. 
Valentini's suggestion on the other hand is to introduce anti-commuting fields as be- 
ables. However, this suggestion meets the problem of identifying a suitable equilibrium 
probability distribution over these anti-commuting fields. Without such a distribution 
it is unclear how to start developing a pilot-wave model. 

One model that has field beables for fermionic fields and that is capable of reproduc- 
ing the standard quantum predictions is the model by Struyve and Westman, which was 
only briefly discussed here. In this model the beables for the fermionic fields are defined 
in terms of the beables for the bosonic degrees of freedom and the wave function. 

On the other hand, taking particle positions as beables seems more succesful for 
fermionic field theories than for bosonic field theories [28-37]. In view of this, it might 
be worthwile to consider a supersymmetric extension of standard quantum field theory. 
In a supersymmetric quantum field theory there is a symmetry between bosons and 
fermions. Maybe this symmetry allows for a more unified approach for introducing 
beables. 
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A On deriving the guidance equation from analogy 
with Hamilton-Jacobi theory 

The general method to find the guidance equation for some configuration q is to consider 
the continuity equation for the quantum mechanical probability density \ip{q,t)\ 2 and 
to define the velocity field for q as v ^ = /\ip\ 2 , where is the probability current 
associated to the density \ip\ 2 . 

For many quantum theories of interest that can be obtained by quantizing a classical 
Lagrangian L(q,q), it is so that the guidance equations q n = j^/l^l 2 are equivalent to 
the equations 

Pn = 7^—, (165) 
dq n 

where p n = dL(q,q)/dq n are the classical momenta and where S is the phase of the 
wave function, i.e. ip = exp(LS) [24]. The latter prescription (165) is inspired by the 
Hamilton-Jacobi formulation of classical mechanics. 

This is for example so for the non-relativistic pilot-wave theory of deBroglie and 
Bohm that was discussed in section 2.1.3. The guidance equations x& = m^V^S* are 
readily obtained by equating the classical momenta, which are given by = m k ± k , with 
'VkS. The pilot-wave theory further yields an equation that resembles the Hamilton- 
Jacobi equation of classical mechanics. However, this resemblance is merely formal. The 
differences are discussed in detail in [17, pp. 131-134]. 

However the prescription (165) to find the guidance equation does not hold in gen- 
eral. We already encountered one example where it doesn't hold, namely in the case 
of the Dirac field, in section 3.3.1. However, the Dirac theory involves both fermionic 
quantization and a singular Lagrangian and it might perhaps be thought that the failure 
of the prescription is due to either of these facts. However, below we consider another 
theory which involves bosonic quantization and a non-singular Lagrangian for which the 
prescription fails too. 

The example we have in mind corresponds to the classical Hamiltonian 



H = ct\p + «2P , (166) 

where ol\ and a<i are strictly positive real constants. The Lagrangian can be obtained 
by performing the Legendre transform, with velocity defined by 

q = ~dp = P + 2a2P ' ^ ^ 



57 



Since d 2 H/dp 2 = 12aip 2 + 2a 2 > 0, the relation above can be inverted to yield the 
momentum p in terms of the velocity q. The Lagrangian is then given by L = qp — H, 
where p is expressed in terms of q. For the subsequent discussion there is no need to 
give the Lagrangian explicitly. 

Hamilton's equation (167) is equivalent to the defining relation of the momentum in 
terms of the Lagrangian, i.e. p = dL/dq. So the above prescription (165) to find the 
guidance equation, yields 



fdS\ 



dS 



4ftl UJ +2a2 dq- 



(168) 



On the other hand, the quantization of the Hamiltonian gives the Schrodinger equa- 
tion: 

for which the corresponding continuity equation reads 



d\^\ 2 + dj^ = Q 



dt 



dq 



(170) 



with 



f = 21m 



4«i 



d 3 ip* ^ dip* d 2 ip\ ^ jj*^ 
1 dq 3 dq dq 2 J dq 



R 2 " 
\dq 



_ 2R d 2 RdS fdR\ 2 dS_ _ R dR^S_ _ l^ 2 d 3 ff 
dq 2 dq \ dq J dq dq dq 2 2 dq 3 



+ 2a 2 R 2 — , 
dq 



(171) 



where we used the polar decomposition ip = Re lS in the last line. The resulting guidance 
equation reads 

. j* 
1 = rrrj 



= 4ai 



+ 2a 2 



dq J dq 2 dq 

dS 



dRydS R _ 1 dRd 2 S _ l d 3 S 
dq J dq ' dq dq 2 2 dq 3 



dq 



(172) 



Hence we see that the correct guidance equation, namely (172), contains extra terms 
not present in (168). 



B Evaluation of the expression (154) 

With 

(")«+= J dVLa\u + \ 2 = ^tt, ( a ) u _=J 



dftalu. 



(173) 
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A M+ « = yVW - = 7^-2. ( 174 ) 



A u _« = ^(« 2 ) u _ - (a>L = y ^ - 2 , (175) 
n the number of particles in V, and rii the number of lattice sites in V, we have 

{A V )^ M = — (n(a) u+ + (m - n)(a) u _) , {A v )m e = («)u_ , (176) 



A* M A V = — A u+ a = — A M _a = A* e A y , (177) 



1 a 1 
— A u+ a = — 

yjn x y/ni 

so that (154) reduces to 



" »,, A -° =V^-2. (178) 



y/ni (a) u+ - (a) u _ tt V 64 

With a the lattice spacing and p the particle number density of the matter under con- 
sideration, and assuming a cubic volume V with edge of length L, we have n = L 3 p and 
ni = (L/a) 3 . Since the right hand side of (178) is of the order one, (178) reduces to 

Lap 213 > 1 . (179) 



References 

[1] L. deBroglie, in Electrons et Photons: Rapports et Discussions du Cinquieme 
Conseil de Physique, Gauthier-Villars, Paris, 105 (1928); English translation: G. 
Bacciagaluppi and A. Valentini, quantum Theory at the Crossroads: Reconsidering 
the 1927 Solvay Conference, Cambridge University Press, Cambridge (2009) and 
arXiv:quant-ph/0609184. 

[2] D. Bohm, Phys. Rev. 85, 166 (1952). 

[3] D. Bohm, Phys. Rev. 85, 180 (1952). 

[4] D. Diirr, S. Goldstein and N. Zanghi, J. Stat. Phys. 67, 843 (1992) and 
arXiv:quant-ph/0308039. 

[5] J.S. Bell, Speakable and unspeakable in quantum mechanics, Cambridge University 
Press, Cambridge (1987). 

[6] J.S. Bell, Epistemological Letters 9, 11 (1976), reprinted in Dialectica 39, 85 
(1985), in J.S. Bell, Speakable and unspeakable in quantum mechanics, Cambridge 
University Press, Cambridge, 52 (1987) and in John S. Bell on The Foundations of 
quantum Mechanics, eds. M. Bell, K. Gottfried and M. Veltman, World Scientific, 
Singapore, 50 (2001). 



59 



[7] J.S. Bell, in Between Science and Technology, eds. A. Sarlemijn and P. Kroes, 
Elsevier, Amsterdam, 97 (1990), reprinted in John S. Bell on The Foundations of 
quantum Mechanics, eds. M. Bell, K. Gottfried and M. Veltman, World Scientific, 
Singapore, 216 (2001). 

[8] V. Allori, S. Goldstein, R. Tumulka and N. Zanghi, Brit. J. Phil. Sci. 59, 353 
(2008) and arXiv:quant-ph/0603027. 

[9] P.N. Kaloyerou, Investigation of the quantum Potential in the Relativistic Domain, 
PhD. Thesis, Birkbeck College, London (1985). 

[10] P.N. Kaloyerou, Phys. Rep. 244, 288 (1994). 

[11] P.N. Kaloyerou, in Bohmian Mechanics and quantum Theory: An Appraisal, eds. 
J.T. Cushing, A. Fine and S. Goldstein, Kluwer, Dordrecht, 155 (1996). 

[12] P.N. Kaloyerou, Found. Phys. Lett. 13, 41 (2000). 

[13] P.N. Kaloyerou, Physica A 355, 297 (2005) and arXiv:quant-ph/0311035. 

[14] P.N. Kaloyerou, J. Phys. A 39, 11541 (2006) and arXiv:quant-ph/0503201. 

[15] D. Bohm, B.J. Hiley and P.N. Kaloyerou, Phys. Rep. 144, 349 (1987). 

[16] D. Bohm and B.J. Hiley, The Undivided Universe, Routledge, New York (1993). 

[17] P.R. Holland, The quantum Theory of Motion, Cambridge University Press, Cam- 
bridge (1993). 

[18] P.R. Holland, Phys. Rep. 224, 95 (1993). 

[19] M.M. Lam and C. Dewdney, Found. Phys. 24, 3 (1994). 

[20] M.M. Lam and C. Dewdney, Found. Phys. 24, 29 (1994). 

[21] T. Takabayasi, Prog. Theor. Phys. 8, 143 (1952). 

[22] P.R. Holland, Phys. Lett. A 128, 9 (1988). 

[23] A. Valentini, On the Pilot-Wave Theory of Classical, quantum and Subquantum 
Physics, PhD. Thesis, International School for Advanced Studies, Trieste (1992), 
online http : //www . sissa . it/ ap/PhD/Theses/valent ini . pdf . 

[24] A. Valentini, in Bohmian Mechanics and quantum Theory: An Appraisal, eds. J.T. 
Cushing, A. Fine and S. Goldstein, Kluwer, Dordrecht, 45 (1996). 

[25] A. Valentini, Hidden variables in modern physics and beyond, in preparation. 



60 



[26] W. Struyve and H. Westman, in quantum Mechanics: Are there quantum Jumps? 
and On the Present Status of quantum Mechanics, eds. A. Bassi, D. Diirr, T. 
Weber and N. Zanghi, AIP Conference Proceedings 844, American Institute of 
Physics, USA, 321 (2006) and arXiv:quant-ph/0602229. 

[27] W. Struyve and H. Westman, Proc. R. Soc. A 463, 3115 (2007) and 
arXiv:0707.3487 [quant-ph]. 

[28] J.S. Bell, CERN preprint CERN-TH. 4035/84 (1984), reprinted in J.S. Bell, Speak- 
able and unspeakable in quantum mechanics, Cambridge University Press, Cam- 
bridge, 173 (1987); Phys. Rep. 137, 49 (1986); in quantum Implications, eds. B.J. 
Hiley and F. David Peat, Routledge, London, p. 227 (1987); in John S. Bell on The 
Foundations of quantum Mechanics, eds. M. Bell, K. Gottfried and M. Veltman, 
World Scientific, Singapore, 159 (2001). 

[29] D. Diirr, S. Goldstein, R. Tumulka and N. Zanghi, J. Phys. A 36, 4143 (2003) 
and arXiv:quant-ph/0208072. 

[30] D. Diirr, S. Goldstein, R. Tumulka and N. Zanghi, Phys. Rev. Lett. 93, 090402 
(2004) and arXiv:quant-ph/0303156. 

[31] D. Diirr, S. Goldstein, R. Tumulka and N. Zanghi, Commun. Math. Phys. 254, 
129 (2005) and arXiv:quant-ph/0303056. 

[32] R. Tumulka and H.-O. Georgii, Interacting Stochastic Systems, eds. J.-D. Deuschel 
and A. Greven, Springer- Verlag, Berlin, 55 (2005) and math.PR/0312326. 

[33] D. Diirr, S. Goldstein, R. Tumulka and N. Zanghi, J. Phys. A 38, Rl (2005) and 
arXiv:quant-ph/0407116. 

[34] S. Colin, arXiv:quant-ph/0301119. 

[35] S. Colin, Phys. Lett. A 317, 349 (2003) and arXiv:quant-ph/03 10055. 

[36] S. Colin, Ann. Fond. Louis de Broglie 29, 273 (2004) and arXiv:quant-ph/0310056. 

[37] S. Colin and W. Struyve, J. Phys. A 40, 7309 (2007) and arXiv:quant-ph/0701085. 

[38] R. Tumulka, J. Phys. A 40, 3245 (2007) and arXiv:quant-ph/0607124. 

[39] S. Twareque Ali and M. Englis, Rev. Math. Phys. 17, 391 (2005) and math- 
ph/0405065. 

[40] P.A.M. Dirac, Lectures On quantum Mechanics, Belfer Graduate School of Science, 
Yeshiva University, New York (1964), reprinted by Dover Publications, New York 
(2001). 

[41] A. Hanson, T. Regge and C. Teitelboim, Constrained Hamiltonian Systems, Ac- 
cademia Nazionale Dei Lincei, Roma (1976). 



61 



[42] K. Sundermeyer, Constrained Dynamics, Lecture Notes in Physics 169, Springer- 
Verlag, Berlin (1982). 

[43] D.M. Gitman and I.V. Tyutin, quantization of Fields with Constraints, Springer- 
Verlag, Berlin Heidelberg (1990). 

[44] M. Henneaux and C. Teitelboim, quantization of Gauge Systems, Princeton Uni- 
versity Press, New Jersey (1991). 

[45] T. Maskawa and H. Nakajima, Prog. Theor. Phys. 56, 1295 (1976). 

[46] C.J. Isham, in Relativity, Groups and Topology II, Les Houches Session XL, eds. 
B.S. Dewitt and R. Stora, North-Holland Physics Publishing, Amsterdam, 1059 
(1984). 

[47] W. Struyve and A. Valentini, J. Phys. A 42, 035301 (2009) and arXiv:0808.0290 
[quant-ph] . 

[48] D. Bohm, Phys. Rev. 89, 458 (1953). 

[49] A. Valentini, Phys. Lett. A 156, 5 (1991). 

[50] D. Diirr and S. Teufel, Bohmian Mechanics, Springer- Verlag, Berlin (2009). 

[51] R. Tumulka, Brit. J. Phil. Set. 58, 355 (2007) and arXiv:quant-ph/0605130. 

[52] A. Valentini and H. Westman, Proc. R. Soc. A 461, 253 (2005) and arXiv:quant- 
ph/0403034. 

[53] S. Colin and W. Struyve, New J. Phys. 12, 043008 (2010) and arXiv:0911.2823 
[quant-ph] . 

[54] D. Bohm and B.J. Hiley, Phys. Rep. 144, 323 (1987). 

[55] D. Diirr, S. Goldstein, R. Tumulka and N. Zanghi, J. Stat. Phys. 116, 959 (2004) 
and arXiv:quant-ph/0308038. 

[56] R. Jackiw, Diverse topics in theoretical and mathematical physics, World Scientific 
Publishing, Singapore (1995). 

[57] B. Hatfield, quantum Field Theory of Point Particles and Strings, Perseus Books, 
Cambridge MA (1991). 

[58] J. Glimm and A. Jaffe, quantum Physics, second edition, Springer- Verlag, New- 
York (1987). 

[59] C.J. Isham, in Recent Aspects of quantum Fields: Proceedings of the XXX Int. 
Universitdtswochen fur Kernphysik, Schladming, Austria, February and March 
1991, eds. H. Mitter and H. Gausterer, Lecture Notes in Physics 396, Springer 
Verlag, Berlin, 123 (1991). 



62 



[60] C.J. Isham, in Integrable Systems, quantum Groups, and quantum Field Theories, 
eds. L.A. Ibort and M.A. Rodriguez, Kluwer Academic Publishers, London, 157 
(1993) and arXiv:gr-qc/9210011. 

[61] P. Colella and O.E. Lanford III, in Constructive quantum Field Theory, eds. G. 
Velo and A. Wightman, Lecture Notes in Physics 25, Springer- Verlag, Berlin, 44 
(1973). 

[62] K. Symanzik, Nucl. Phys. B 190 [FS3], 1 (1981). 

[63] M. Liischer, Nucl. Phys. B 254, 52 (1985). 

[64] F. Cooper and E. Mottola, Phys. Rev. D 36, 3114 (1987). 

[65] S.-Y. Pi and M. Samiullah, Phys. Rev. D 36, 3128 (1987). 

[66] M. Liischer, R. Narayanan, P. Weisz and U. Wolff, Nucl. Phys. 53844, 168 (1992). 
[67] J. Kogut and L. Susskind, Phys. Rev. D 11, 395 (1975). 
[68] J.B. Kogut, Rev. Mod. Phys. 51, 659 (1979). 

[69] S. Saunders, in From Physics to Philosophy, eds. J. Butterfield and C. Pagonis, 
Cambridge University Press, Cambridge, 71 (1999). 

[70] K. Berndl, D. Diirr, S. Goldstein and N. Zanghi, Phys. Rev. A 53, 2062 (1996) 
and arXiv:quant-ph/9510027. 

[71] D. Diirr, S. Goldstein, K. Miinch-Berndl and N. Zanghi, Phys. Rev. A 60, 2729 
(1999) and arXiv:quant-ph/9801070. 

[72] A. Valentini, Phys. Lett. A 228, 215 (1997). 

[73] L.A. Gergely, Ann. Phys. 298, 394 (2002). 

[74] K. Baumann, II Nuovo Cimento B 96, 21 (1986). 

[75] R. Weingard, in Perspectives on quantum Reality, ed. R. Clifton, Kluwer Academic 
Publishers, Dordrecht, 195 (1995). 

[76] N. Huggett and R. Weingard, in Conceptual issues in quantum field theory, ed. 
T.Y. Cao, Cambridge University Press, Cambridge, 287 (1999). 

[77] R.J. Glauber, Phys. Rev. 131, 2766 (1963). 

[78] I. Schmelzer, Found. Phys. 40, 289 (2010) and arXiv:0904.0764 [quant-ph]. 

[79] W. Struyve, "Spontaneous symmetry breaking and the Higgs mechanism: Lifting 
the veil of gauge symmetry" , in preparation. 



63 



[80] F.J. Belinfante, A survey of hidden-variable theories, Pergamon Press, Oxford 
(1973). 

[81] P.A.M. Dirac, The Principles of quantum Mechanics, fourth revised edition, 
Clarendon Press, Oxford (1967). 

[82] V.N. Gribov, Nucl. Phys. B 139, 1 (1978). 

[83] I.M. Singer, Comm. Math. Phys. 60, 7 (1978). 

[84] A. Chodos and V. Moncrief, J. Math. Phys. 21, 364 (1980). 

[85] J. A. Wheeler, in Relativity, Groups, and Topology, eds. B. DeWitt and C. DeWitt, 
Gordon and Breach Science Publishers, New York, 315 (1964). 

[86] J. A. Wheeler, in Battelle Rencontres, eds. CM. DeWitt and J. A. Wheeler, W. A. 
Benjamin, Inc., New York, 242 (1968). 

[87] R. Healey, Gauging What's Real, Oxford University Press, Oxford (2007). 

[88] E.S. Abers and B.W. Lee, Phys. Rep. 9, 1 (1973). 

[89] D. Karabali and V.P Nair, Nucl. Phys. B 464, 135 (1996) and arXiv:hep- 
th/9510157. 

[90] PR. Holland, Phys. Rep. 169, 293 (1988). 

[91] J.S. Bell, Rev. Mod. Phys. 38, 447 (1966), reprinted in J.S. Bell, Speakable and 
unspeakable in quantum mechanics, Cambridge University Press, Cambridge, 1 
(1987) and in John S. Bell on The Foundations of quantum Mechanics, eds. M. 
Bell, K. Gottfried and M. Veltman, World Scientific, Singapore, 1 (2001). 

[92] J.S. Bell, in Foundations of quantum Mechanics. Proceedings of the International 
School of Physics 'Enrico Fermi', course IL, Academic, New York, 171 (1971), 
reprinted in J.S. Bell, Speakable and unspeakable in quantum mechanics, Cam- 
bridge University Press, Cambridge, 29 (1987) and in John S. Bell on The Foun- 
dations of quantum Mechanics, eds. M. Bell, K. Gottfried and M. Veltman, World 
Scientific, Singapore, 22 (2001). 

[93] P. Jordan and E. Wigner, Z. Phys. 47, 631 (1928). 

[94] R. Floreanini and R. Jackiw, Phys. Rev. D 37, 2206 (1988). 

[95] C. Kiefer and A. Wipf, Ann. Phys. 236, 241 (1994) and arXiv:hep-th/9306161. 

[96] F.A. Berezin, The Method Of Second quantization, Academic Press, New York 
(1966). 



64 



[97] F.A. Berezin, Introduction to Superanalysis, D. Reidel Publishing Company, Dor- 
drecht (1987). 

[98] B. DeWitt, Supermanifolds, second edition, Cambridge University Press, Cam- 
bridge (1992) . 

[99] J. Hallin and P. Liljenberg, Phys. Rev. D 52, 1150 (1995) and arXiv:hep- 
th/9412188. 

[100] W. Struyve, The de Broglie-Bohm pilot-wave interpretation of quantum theory, 
PhD. Thesis, Ghent University, Ghent (2004) and arXiv:quant-ph/0506243. 

[101] J.L. Martin, Proc. R. Soc. A 251, 536 (1959). 

[102] F.A. Berezin and M.S. Marinov, Ann. Phys. 104, 336 (1977), reprinted in F.A. 
Berezin, Introduction to Superanalysis, D. Reidel Publishing Company, Dordrecht 
(1987). 

[103] R. Casalbuoni, Nuovo Cimento 33, 115 (1976). 
[104] R. Casalbuoni, Nuovo Cimento 33, 389 (1976). 

[105] A. Bassi and G. Ghirardi, Phys. Rep. 379, 257 (2003) and arXiv:quant- 
ph/0302164. 

[106] R. Gambini and J. Pullin, Loops, Knots, Gauge Theories and quantum Gravity, 
Cambridge University Press, Cambridge (1996). 

[107] C. Rovelli, quantum Gravity, Cambridge University Press, Cambridge (2004). 

[108] T. Thiemann, Modern Canonical Quantum General Relativity, Cambridge Uni- 
versity Press, Cambridge (2007). 



65 



